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[ Abstract. Between 1996 and 2002, we have carried out a spectral monitoring program for the Seyfert galaxy 

NGC 5548 with the 6 m and 1 m telescopes of SAO (Russia) and with the 2.1 m telescope of Guillermo Haro 
' Observatory (GHO) at Cananea, Mexico. High quality spectra with S/N> 50 in the continuum near Ha and H/3 

were obtained, covering the spectral range ~(4000 - 7500) A with a (4.5 to 15) A-resolution. We found that both 
the flux in the lines and the continuum gradually decreased, reaching minimum values during May- June 2002. 
In the minimum state, the wings of H/3 and Ha became extremely weak, corresponding to a Syl.8 type, not to 
a Syl, as observed previously when the nucleus was brighter. The line profiles were decomposed into variable 
and constant components. The variable broad component is well correlated with the continuum variation. It 
consists of a double peaked structure with radial velocities ~ ±1000 km/s relative to the narrow component. A 
constant component, whose presence is independent of the continuum flux variations, shows only narrow emission 
lines. The mean, rms, and the averaged over years, observed and difference line profiles of H/3 and Ha reveal the 
same double peaked structure. The relative intensity of these peaks changes with time. During 1996, the red peak 
was the brightest, while in 1998 - 2002, the blue peak became the brighter one. Their radial velocities vary in 
the ~ (500 - 1200) km/s range. In 2000 - 2002 a distinct third peak appeared in the red wing of Ha and H/3 
line profiles. The radial velocity of this feature decreased between 2000 and 2002: from the observed profiles, 
from ~ +(2500 - 2600) km/s to ~ +2000 km/s and is clearly seen on the difference profiles. The fluxes of the 
various parts of the line profiles are well correlated with each other and also with the continuum flux. The blue 
and red parts of the line profiles at the same radial velocities vary in an almost identical manner. Shape changes 
of the different parts of the broad line are not correlated with continuum variations and, apparently, are not 
related to reverberation effects. Changes of the integral Balmer decrement are, on average, anticorrelated with 
the continuum flux variations. This is probably due to an increasing role of collisional excitation as the ionizing 
flux decreases. The behavior of the Balmer decrement of the various parts of the line profiles was different in 
1996 - 2000 as compared with the 2001 behavior. Our results favor the formation of the broad Balmer lines in 
a turbulent accretion disc with large and moving "optically thick" inhomogeneities, capable of reprocessing the 
central source continuum. 
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1. Introduction 

An important question in the study of active galactic nu- 

clei (AGN) is the nature of the "central engine" . A popular 

Send offprint requests to: assumption is that the nuclear activity is caused by accre- 

A.I. Shapovalova e-mail: ashap@sao.ru tion of gas onto a supermassive black hole (Rees 111)841 
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Begclman ll985fl . The basic energy release of an AGN oc- 
curs very close to the nucleus (r < 0.001 pc), as an UV and 
X-ray continuum most probably produced by a geometri- 
cally thin accretion disc. Then, broad emission lines are 
produced in a zone (the BLR) that reprocesses a fraction 
of the central UV-X continuum. A zone located further 
out (r > 0.001 pc). The BLR is filled with gas obviously 
linked with the accretion process. It is therefore important 
to know its structure and kinematics, in order to gain in- 
sight into the central engine. However, even for nearby 
objects, the typical angular size of the BLR corresponds 
to < 0.001 arcsec, hence we will have to wait for the avail- 
ability of more sensitive optical interferometers to resolve 
it. Fortunately, another methods exist to study the BLR 
structure. 

It is well known that AGNs vary in luminosity on 
time scales from years to hours, over the entire wavelength 
range from radio to X-ray or 7-ray. In particular, the flux 
in the broad emission lines varies in response to changes 
in the ionizing continuum with short time delays (days to 
weeks for Seyfert galaxies), due to light-travel time effects 
within the BLR. If the BLR gas has systematic motions 
such as infalling, outflowing, circular motions, etc., then 
the profiles of the broad emission lines must vary in a way 
related with the geometry and the kinematics of the gas in 
this region, and with the processes of gas relaxation that 
follow the changes in the ionizing flux (Bahcall, Kozlovsky 
& Salpeter 1972; Bochkarev & Antokhin Blandford 
& McKeeHHES Antokhin & Bochkarev HMty . 

Studying the correlations between the flux changes 
in the continuum and in the broad emission line pro- 
files, one can obtain a "map" of the geometrical and dy- 
namical structure of the BLR. This method is known 
as "Reverberation mapping" (see Peterson 119931 and ref- 
erences therein). Important progress in the understand- 
ing of the BLRs was achieved as a result of multiwave- 
length monitoring campaigns within the framework of the 
"International AGN Watch" , a consortium organized to 
study several Seyfert galaxies (Peterson et al. 11999(1 . A 
large amount of data has been obtained in the multiwave- 
length monitoring of the Seyfert 1 galaxy NGC 5548, in- 
cluding its continuous monitoring in the optical range for 
13 years (1988 - 2001) (Clavel et al. ITMT1 Peterson et 
al. HTM! IT39H IT5391 1271H21 Dietrich et al. ITM31 Korista 
et al. H995J) . These investigations have given the following 
results: 

1. The response time of the H/3 line to continuum vari- 
ations varies from year to year (from ~ 8 to 26 days) 
and it is correlated with the average continuum flux. 

2. The high ionization lines have shorter response times 
to continuum variations than the low ionization ones, 
indicating the presence of ionization stratification 
along the radius of the broad line region. 

3. The optical and ultraviolet continua vary almost simul- 
taneously (i.e. without showing a significant delay). 

4. The UV/optical continuum becomes "harder" as it 
gets brighter. 



The results of the study of the radial velocity field in the 
BLR of NGC 5548 are still ambiguous. The analysis of 
the "AGN Watch" spectra reveals a complex behaviour. 
Crenshaw and Blackwcll (1990) found that the red wing 
of C iv A 1550 responds faster to changes of the continuum 
than the blue wing does, implying important radial mo- 
tions (infalling). Later on, Korista l|1994fl observed that 
both wings presented the same delay with respect to con- 
tinuum light variations, implying Keplerian rotation or 
turbulent symmetrical motions. The absence of important 
radial motions in the BLR was further confirmed by the 
analysis of the transfer function (TF) of the H/3 (Wanders 
& Peterson UMB) and CivA1550 (Wanders et al. IT§55|> 
emission lines. It was also found that the broad wings 
of the emission lines respond faster to continuum varia- 
tions than the line cores (Clavel et al. 119911 KoristaQ r 994 
Korista et al. 119951 Kollatschny & Dietrich 1996, Wanders 
& Peterson [1996). In several studies, the presence of a 
multi-component structure of the BLR with distinct phys- 
ical characteristics is required to account for the observed 
broad line profile variability in this object (Peterson 1987 
Stirpe et al. HUggl Stirpe & de Bruyn ITMT1 Sergeev et 
al. 119941 Wanders and Peterson ll996|l . Consequentely, the 
observed profiles in NGC 5548 have been interpreted in 
the framework of different models, such as: binary black 
holes (Peterson et al. I1987[) . accretion disc (Stirpe et 
al. EES Rokaki & Boisson lT55§jl . BLR from the clouds 
rotating about the massive central gravitational source in 
occasionally inclined Keplerian orbits and illuminated by 
an anisotropic continuum source (Wanders et al. [1995 
Wanders & Peterson 1 1996|l . 

Further study of the broad emission line profile changes 
on longer time scales, may allow us to prove or disprove 
some of the previously advanced hypothetical scenarios. 
This is the purpose of the present paper. 

In this paper we present the results of an optical spec- 
tral study of NGC 5548 for the 1996 - 2001 period, includ- 
ing part of our 2002 data (see section |2~T|) . Some partial 
results of our monitoring campaign were reported earlier 
(Shapovalova et al. I2001al l2MTbl [MUt- In section El we 
discuss the observations and data processing. In section[21 
we present the analysis of the H/3 and Ha line profile vari- 
ability. The decomposition of the profiles into constant 
and variable components, the mean and rms spectra, are 
discussed. The behaviour of the radial velocities of spec- 
tral features in the broad line profiles, is studied both from 
the observed profiles and from the difference profiles. We 
use words "bump" for wide features and "peak" for their 
tops (maxima). The correlation analysis between fluxes 
and shapes of the different parts of the line profiles is pre- 
sented. The behaviour of the Balmer decrement is studied. 
In Section^] we summarize our results and compare them 
with those of other researchers. Possible interpretations 
are discussed in Section and conclusions are listed in 
Section In Appendix ^ a modified spectrum scaling 
method adopted in this paper is presented. 

Note that all these data are also included in the AGN 
Watch database and are publicly available. 
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Table 1. Sources of Spectroscopic Observations 



Source 


Code 


Tel. and Equip. 


Aperture 


1 


2 


3 


4 


SAO (Russia) 


LI 


lm+UAGS 


4.2" x 19.8" 


SAO (Russia) 


LI 


lm+UAGS 


8.0" x 19.8" 


SAO (Russia) 


L 


6m+UAGS 


2.0" x 6.0" 


Guillermo Haro 


GH 


2.1m+B&C 


2.5" x 6.0" 


Obs. (Mexico) 









2. Observations and data reduction 

2.1. Optical Spectroscopy 

We report the spectral observations of NGC 5548 carried 
out between 1996 Jan 14 (Julian date= JD 2 450 097) and 
2001 Aug 9 (JD 2 452131) during 113 nights. Our anal- 
ysis is based on those spectra. However for the studies 
described in Sections 13.11 13.2.31 and 13.2.41 (for studing the 
behavior of the peaks and bumps) we have also used spec- 
tra taken on June 4, 2002, and May 15, 17, 2002, when 
NGC 5548 was in a minimum activity state and the 2002 
annual averages of both the observed and difference pro- 
files for Ha and H/3. 

Optical spectra of NGC 5548 were obtained with the 
6 m and 1 m telescopes of SAO (Russia, 1996 - 2002) 
and at INAOE's 2.1 m telescope at the Guillermo Haro 
Observatory (GHO) at Cananea, Sonora, Mexico (1998 - 
2002). These were obtained with long slit spectrographs 
equipped with CCDs. The typical wavelength range cov- 
ered was from 4000 A to 7500 A, the spectral resolution 
was 4.5 - 15 A, and the S/N ratio was > 50 in the contin- 
uum near Ha and H/3. Spectrophotometric standard stars 
were observed every night. The informations on the source 
of spectroscopic observations are given in Tabled 1 - the 
source (Observatory); 2 - a code assigned to each tele- 
scope+cquipmcnt, used throughout this paper (the code 
was chosen in accordance with the monitoring campaigns 
of NGC 5548 (Dietrich et al. 120011 Peterson et al. 12002)) : 
3 - the telescope aperture and the spectrograph; 4 - the 
projected spectrograph entrance apertures (the first di- 
mension is the slit-width, and the second one is the slit- 
length) . 

Table 2. Log of the spectroscopic observations 

This table is only available in electronic form at the CDS. It 
contains the following information for 116 dates: Columns: 1 
- number; 2 - UT date; 3 - Julian date; 4- code according to 
Table 1; 5 - projected spectrograph entrance apertures; 6 - 
wavelength range covered; 7 - spectral resolution; 8 - mean 
seeing; 9 - position angle (PA) in degrees; 10 - signal to noise 
ratio in the continuum (5160 - 5220) AA near H/3 and (6940 - 
7040) AA near Ha. 

The spectrophotometric data reduction was carried 
out either with the software developed at the SAO RAS 
by Vlasyuk ( 1993), or with IRAF for the spectra obtained 



in Mexico. The image reduction process included bias 
subtraction, flat-field corrections, cosmic ray removal, 2D 
wavelength linearization, sky spectrum subtraction, stack- 
ing of the spectra for every set-up, and flux calibration 
based on standard star observations. 

2.2. Absolute calibration of the spectra 

Even under good photometric conditions, the accuracy 
of spectrophotometric measurements is rarely better than 
10%. Thus the standard technique of flux calibration, by 
means of comparison with stars of known spectral energy 
distribution, is not good enough for the study of AGN vari- 
ability. Instead, we use the fluxes of the narrow emission 
lines which are known to be non- variable on time scales of 
tens of years in most AGN. Consequently, the bright nar- 
row emission lines can be adopted as internal calibrators 
for scaling AGN spectra (Peterson 1993). So, we assume 
that the flux of the [O in] A 5007 line remains constant dur- 
ing the interval covered by our observations. All blue spec- 
tra of NGC 5548 are scaled to a constant flux value of 
.F([Oiii]A5007) = 5.58 x 10" 13 ergs" 1 cm~ 2 determined 
by Peterson et al. I|1991[l and corrected for aperture ef- 
fects as described below. The scaling of the blue spectra 
was carried out using a variation on the method of Van 
Groningen & Wanders H1992|l . described in AppendixlAl 
This method allows to obtain a homogeneous set of spec- 
tra with the same wavelength calibration and the same 
[OIII] A5007 flux value. 

The spectra obtained with the GHO 2.1 m telescope 
(Mexico) with a resolution of 15A. They contain both the 
Ha and H/3 regions, and were scaled using the [O in] A 5007 
line. 

Most of spectra from the 1 m and 6 m SAO telescopes 
were obtained separately in the blue (H/3) and red (Ha) 
wavelength intervals, with a resolution of 8 - 9 A. Usually, 
the red edge of the blue spectra and the blue edge of the 
red spectra overlap in an interval of ~ 300 A. Therefore, 
the most of the red spectra were scaled using the over- 
lapping continuum region with the blue ones, which were 
scaled with the [O in] line. In these cases the scaling un- 
certainty is about 5%. 

However, for 10 red spectra the scaling of the contin- 
uum by this method was not possible, this due to several 
reasons: ie. some spectra were obtained with a higher res- 
olution (~ 5A) and did not overlap with the blue spectra; 
or the some spectrum ends were distorted by the reduction 
procedure; or blue spectra were not taken on that night. 
These 10 spectra were scaled using the integral flux in the 
narrow emission lines in the Ha, region: [Nn]AA 6548, 
6584 and [Sn]AA 6717, 6731. To this purpose, on the red 
spectra we located a linear continuum through points that 
are free from the absorption lines (6120 A and 7020 A) in 
20 Awindows. 

After continuum subtraction, we obtained the best 
gaussian fit to the Ha profile through a blend of 
7 emission components: very broad Ha (FWHM~ 
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Table 3. Fluxes of Ha Gaussian components 



6500 6600 6700 6600 6900 

Wavelength (A) 

Fig. 1. Gaussian deconvolution of the Ha blend in sev- 
eral components; thick solid line: the observed profile; thin 
solid line: the narrow component of Ha, [Nn]AA6548, 
6584, [Sn]AA6717, 6731; dashed line: very broad, blue 
broad, and red broad components of Ha; points: the resid- 
ual of the observed profile and the sum of all gaussian 
components. Our analysis do not include a small hump in 
the left side ([Oi]A6364). 



10000 km/s), blue broad Ha, red broad Ha, narrow 
Ha, [N n]AA 6548+6583 (double gaussian function with 
I(6584)/I(6548)=3), [Sn]A6717 and [Sn]A6731. Fig. □ 
shows an example of a Gaussian fit to the Ha blend. 

Table [3] lists our mean values of the flux in the narrow 
components derived from the scaled red spectra, jointly 
with the result obtained by Dietrich et al. ill 003 ■ from the 
spectra taken in 1088 - 1989. These authors obtained a 
good Gaussian fit to Ha with only 5 components: namely, 
very broad Ha, broad Ha, narrow Ha, [N n] AA 6548, 6584. 
With our spectra we tried to carry out a similar fitting to 
the same 5 components, but we ended up with large resid- 
uals (~ 10 - 20%) in the region where the blue and red 
bumps appear. Apparently, on the spectra in Dietrich et 
al. I|1993f) . such features were not as bright as in our spec- 
tra. However, in any profile fit, the narrow components 
should remain constant. 

We can see on Table |3] that the fluxes obtained for the 
narrow lines are consistent within the errors with the val- 

EH 



ues obtained by Dietrich et al. (jX993|. Using the integral 
fluxes for the narrow lines from Table we scaled the 10 
red spectra. 



Component 


Flux* 






mean (GHO+SAO) 


Dietrich et al. 






(1993) 


Ho^narrow 


0.350 ± 0.040 


0.344 ± 0.045 


[Nil] A A 6548, 6584 


0.160 ±0.036 


0.182 ± 0.022 


[Sn]AA6717, 6731 


0.097 ±0.021 




Sum (H Qn , [Nil]) 


0.510 ±0.054 


0.526 ± 0.050 


Sum narr. all 


0.607 ± 0.064 





|! F([Oiii]A5007)=1.0 



2.3. Intercalibration of the spectral data 

^From the scaled spectra we determined the average flux 
in the continuum at the observed wavelength ~ 5190A (or 
at - 5100A in the rest frame of NGC 5548, z=0.0167), by 
means of flux averages in the bandpass (5180 - 5200) A. 
For the determination of the observed fluxes of H/3 and 
Ha, it is necessary to subtract the continuum. To this 
goal, a linear continuum was located through windows of 
20 A located at 4790 A and 5 170 A for the H/3 region, 
and at 6120 A and 7020 A for the Ha region. After the 
continuum subtraction, we defined the observed fluxes in 
the lines in the following wavelength intervals: (4795 - 
5018) A for H/3 (the interval is similar to that in Peterson 
et al. VM^ . and (6500 - 6800) A for Ha (the interval is 
like that in Dietrich et al. I200ip . 

All fluxes were corrected for aperture effects because, 
while the BLR and non-stellar continuum are effectively 
point-like sources (<C 1"), the NLR is an extended one 
(> 2"). Consequently, the measured NLR flux depends 
on the size of the spectrograph's entrance aperture (see 
Peterson et al. 119951 for a detailed discussion) . In order 
to correct our fluxes for aperture effects, we determined a 
point-source correction factor (</?) given by: 



F(H/3) = ip • F A5 oo7 



( F(H/3) 
VF([Oiii]A5007) 



(1) 



where FA5007 is the absolute flux in the [O ill] A 5007 line, 
and the value in brackets is the H/3 to [O in] A 5007 ob- 
served flux ratio. 

The light contribution of the host galaxy to the con- 
tinuum depends also on the aperture size. The continuum 
fluxes F\(5190) were corrected for different amounts of 
host-galaxy contamination, according to the following ex- 
pression (see Peterson et al. I1995fl : 



F A (5190A) =ip-F 



A5007 



/ F A (5190A) \ 
VF([Om]A5007) J 



-G, (2) 



obs 



where F A 5oo7 is the absolute flux in the [Om]A5007 line 
and the value in brackets is the continuum to [O in] A 5007 
observed flux ratio, G being an aperture dependent cor- 
rection factor to account for the host galaxy light. The 
case LI (Tabled, which corresponds to a relatively large 
aperture (4.2" x 19.8"), was adopted as a standard (i.e. 
ip = 1.0, G=0 by definition). The corrections (p and G 
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Table 4. Flux Scale Factors for Optical Spectra 



Sample 


Aperture 


Point-Source 


Extended Source 






Scale factor 


Correction 






H> 


G* 


LI 


4.2" x 19.8" 


1.000 


0.000 


LI 


8.0" x 19.8" 


1.089 ± 0.031 


0.84 ±0.45 


GH 


2.5" x 6.0" 


1.009 ± 0.021 


-2.76 ±0.11 


L 


2.0" x 6.0" 


1.036 ± 0.048 


3.85 ±0.86 



* in units of 10 15 erg s 1 cm 2 A 



were defined for every aperture via the comparison of a 
pair of observations separated in time by 1-3 days. This 
means that the real variability on shorter times (< 3 days) 
will be suppressed by the procedure of data recalibration. 
Point-source correction factors ip and G values for various 
samples are given in Tabled Using these factors, we recal- 
ibrate the observed fluxes in the lines and continuum to a 
common scale corresponding to an aperture of 4.2" x 19.8". 

Table 5. Observed continuum, H/3 and Ha fluxes. 

Table |^| where we list our results of fluxes measurement, is 
only available in electronic form at the CDS and contains the 
following information. Columns: 1 - UT-Date; 2 - Julian date; 
3- a telescope code, according to Tabled 4 - F(cont), the 
continuum flux at 5190 A (in units of 10 -15 ergs -1 cm -2 A -1 ), 
reduced to the 1 m telescope aperture 4.2" x 19.8"; 5 - e c , the 
estimated continuum flux error; 6 - F(H./3), the H/3 total flux 
(in units of 10 -13 ergs -1 cm -2 ); 7 - £h/3, the H/3 flux error; 8- 
F(Ra), the Ha total flux (in units of 10 -13 ergs -1 cm -2 ); 9 - 
£Ha, the Ha flux error. 

The fluxes listed in Table were not corrected for the 
contributions of the narrow-line emission components of 
H/3, Ha, and [Nn]AA6548, 6584. These are constant and 
should not influence a broad line variability study. The 
mean error (uncertainty) in our flux determination for 
both, the H/3 and the continuum, is ~ 3%, while it is 
~ 5% for Ha. These quantities were estimated by com- 
paring our results from the spectra obtained within time 
intervals shorter than 3 days. 

2.4. Subtraction of the narrow emission line 
contribution 

In order to study the broad components of hydrogen lines 
showing the main BLR characteristics, one must remove 
the narrow component of these lines and the forbidden 
lines from the spectra. To this purpose, we constructed 
spectral templates for the H/3 and the Ha blends using 
a Gaussian fit to the higher spectral resolution (~ 8 A) 
profiles observed near the minimum light state (for de- 
tails see section 12.21 and Fig. ^) . The obtained template 
spectra are the sum of the following gaussian compo- 
nents: for H/3, the narrow component of H/3, [O ni]AA 4959, 
5007; for Ha, the narrow component of Ha, [Nn]AA6548, 



6584 and [Sn]AA6717, 6731. The flux values obtained 
for the narrow components of H/3 and Ha in the tem- 
plate spectra are: F(H/3)„ = (0.131 ± 0.012)F(A5007); 
F(Ha)„ = (0.398 ±0.02).F(A5007), respectively. Then, we 
scaled the blue and red spectra according to our scaling 
scheme (see Appendix^), taking the template spectrum 
as a reference. Our template spectrum and any individ- 
ual observed spectrum are thus matched in wavelength, 
reduced to the same resolution, and then subtracted from 
one another. After subtraction of the narrow components, 
the spectra of the Ha and H/3 broad lines is reduced to the 
aperture 4.2" x 19.8" (equation QJ, using ip values listed 
in Table H (section |2~3|) . 

Since the scaling is done from the line [O in] A 5007, we 
had to use the FWHMs obtained from the corresponding 
blue spectrum to reduce the Ha region. But the FWHM 
of [Nn]AA6548, 6584 and [Sn]AA6717, 6731 lines is usu- 
ally somewhat smaller (by ~ 5 - 10%) than the FWHM 
of [O in] A 5007. Therefore, the subtraction of these compo- 
nents is rather poor when one has large spectral resolution 
differences between the template and individual spectra. 

Another way to remove the narrow emission lines, con- 
sists in subtracting the spectra from each other, i.e. ob- 
taining difference spectra. The H/3 and Ha difference pro- 
files are thus obtained by subtracting the minimum ac- 
tivity state spectrum from the individual spectra. For a 
good subtraction, it is necessary for the spectra to have 
similar spectral resolution. These questions are solved via 
the method mention in Appendix [S] where a spectrum 
in minimum state is used as a reference spectrum. The 
difference spectra will be discussed hereinafter. 

3. Data analysis 

3.1. Variability of Hj3 and Ha emission lines and of 
the optical continuum 

In Table[S]it is apparent that both, the continuum and per- 
mitted line emission fluxes, decreased continuously from 
maximum values in 1998 to minimum ones in 2002. The 
maximum amplitude ratios of the flux variations during 
this period were: for H/3 line - ~ 4.7; for Ha - ~ 3.4; and 
for the AA5190 A continuum - ~ 2.5. Photometric data by 
Doroshenko et al.( l2001|l and Spiridonova O.I.( :20Q2j) also 
yield a maximum amplitude of the flux variations ratio 
~ 2.5 in the V band in this period. Thus, the variations 
inferred from the spectral and broad band photometry 
of the continuum are in excellent agreement. This fact is 
in turn indicative of a constant flux in the line [Oin]A 
5007, a key assumption in our spectral scaling scheme. It 
is worth mentioning, that the galactic bulge light contri- 
bution to our spectra was not subtracted. Therefore, the 
derived maximum amplitude ratio for the continuum flux 
changes is, consequently, smaller than the ones derived 
for emission lines. In 2002, the flux in the lines and con- 
tinuum reached a minimum value by mid May or early 
June. In fig.|21we plot spectra of the high and low activity 
states, obtained in June 26 1998, and June 4 2002, respec- 
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tivcly. There one can see that in the low activity state, the 
flux in the continuum decreased by a large factor (~ 2.5 
times), while the slope of the continuum in the blue, be- 
came significantly flatter, showing a spectral index > 2 to 
be compared with the ~ 1.0 a value observed at the high 
activity state. Also, at minimum activity state, the emis- 
sion wings of H/3 and Ha became extremely weak. These 
profiles correspond to a Syl.8 type and not to a Syl, as 
observed in maximum light (i.e. the spectral type of the 
object had suffered a dramatic change!). 

In Fig. |2| we present the light curves of the H/3 and 
the Ha emission lines and the combined continuum. The 
flux in the lines and spectral continuum are those listed 
in Table |SJ The combined continuum data includes both 
the A5190 spectral data and some V band photometry 
from Doroshenko et al. 1|2()01|) and Spiridonova (2002). 
The V band data were converted into A5190 continuum 
flux values through the expression given by Dietrich et al. 
(20011: 



(3) 



logF A (5500) = -0.4 my - 8.439. 



A comparison of the spectral continuum at A5190A, 
-F C (A5190) with some simultaneous observations made in 
the V band leads to the following transformation equation: 



F C (A5190) = (1.120 ±0.064) x F X (V) 
-(0.608 ±0.776), 



(4) 



(correlation coefficient r= 0.95). Here the fluxes are given 
in units of 10~ 15 erg cm _2 s _1 A -1 . One sees on Fig. 3 
that the H/3 and Ha fluxes vary approximately in the 
same way. We have compared our light curves of H-beta 
and of the continuum with those obtained by Peterson 
et al. H2002|l from the AGN- Watch program. The com- 
parison shows that our H-beta fluxes coincide within 
the uncertainties with the AGNW data. The continuum 
fluxes differ because of the different apertures, which are 
4.2" x 19.8" for our data and 5.0" x 7.5" for the AGNW 
data. The correction for the different apertures necessary 
to convert our data to Peterson et al.'s data l|2002(l is 
F(cnt)pct.=F(cnt)our -2.5 (the fluxes are in units 10 -15 
erg cm _2 s _1 A -1 ). With the introduction of such a cor- 
rection, our continuum fluxes coincide with the AGNW 
data within the uncertainties. 



3.2. Ha and Hf3 profile variations 
3.2.1. H/3 and Ha profile decomposition 



As shown by Sergeev et al. Q1994JI . if a line profile is com- 
posed of only a variable and a constant part, it is pos- 
sible to separate these components. The constant com- 
ponent is present in the form of narrow emission lines, 
which are formed far away from the nucleus and are in- 
dependent of changes in the continuum emitted by the 
central source. The variable component is formed in the 
BLR and is strongly dependent on the flux of the cen- 
tral source through reverberation. We have applied the 



1998 Jun26; JD2450991 (high state) 
.2002 Jun4; JD2452430 (low state) 
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Fig. 2. The spectra of NGC 5548 corresponding to the 
high activity state (top) and to the low activity state 
(bottom). The absissae give the observed wavelengths 
(z=0.0167). The ordinates give the fluxes in units of 

-2 s- 1 A" 1 . 



10 15 erg cm 



1 996 



99J 



2000 




500 1000 1500 2000 
JD-2450000, clays 

Fig. 3. The light curves of H/3, Ha, and the contin- 
uum at the observed wavelength A5190A. The nar- 
row components were not subtracted. Fluxes in lines 

and in continuum in 
. The continuum is the "combined 
continuum" obtained from the spectra and photometric 
points in the V band in 1996 - 2001 



are in units of 10 13 erg cm 2s 1 
10~ 15 erg cm~ 2 s _1 A 



two-component profile decomposition method to spectra 
obtained in 1996 - 2001. We assume that the fluxes in 
the variable part of the emission lines are correlated lin- 
early with the continuum flux, as shown by Peterson et al. 
(2002] Fig. 3). An underlying continuum was subtracted 
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from each of the scaled spectra. The region occupied by 
an emission line is divided into equidistant narrow spec- 
tral intervals (5 A bins). For every bin a light curve is con- 
structed. A linear regression for each of these light curves 
and that of the continuum was computed. This is done by 
taking into account the mean delay time between them 
(~ 20 days, Peterson et al. ( 20021)) ■ The mtroduct ion of 
the more precise time delays obtained for each year by 
Peterson et al. (2002) for H/3 does not change the profile 
shape of the variable component, within the errors marked 
in Fig.0] Since the exact year-averaged time delay for Ha 
is not known, we used the same average time delay for the 
both lines. The variable part in the profile is computed as 
the increment of line flux per unit flux increase of the con- 
tinuum. The non-variable part of the profile is estimated 
by the extrapolation of the line flux in every bin to a zero 
continuum flux value. 

This scheme allows us to separate the part of the line 
profiles that is linearly related to the continuum light 
changes from those which remain constant in time. It also 
provides a scheme to estimate the H/3 and Ha narrow 
line fluxes in an independent way from that described in 
Sect. 12.41 In any case both methods yielded very similar 
results. 

Fig-Elshows the variable (top) and constant (bottom) 
components of the H/3 (left) and Ha (right) lines. The cor- 
relation coefficient of the variable component with the con- 
tinuum across the emission lines is plotted in the central 
panels. Both variable components present a double-peaked 
structure with maxima at radial velocities ~ ±1000 km/s. 
The variable part of the H/3 and Ha profiles between 
~ —4000 km/s and ~ +5000 km/s show a highly corre- 
lated (r~ 0.8) response to changes in the continuum. The 
constant component panels contain mainly the narrow line 
emission of H/3, [Om] A 4959, 5007 A (left) and of Ha, 
[Nil] AA6548, 6584, [Sn] AA6717, 6731 and [Oi] A 6364 
(right). 

3.2.2. Mean and Root-Mean-Square Spectra 

The comparison between an average and root-mean- 
square (rms) spectrum provides us a good measure of the 
profile variability. Average and rms H/3 and Ha profiles 
were calculated after removing the continuum and the nar- 
row lines from the profiles (see Sections 12.31 and l2.4fl . 

The mean H/3 and Ha profiles and the absolute rms 
variations per unit wavelength are shown in Figs. It is 
clear that both profiles present a double-peaked structure 
with maxima at radial velocities ~ ±1000 km/s relative 
to the narrow components. On the mean and rms Ha and 
H/3 profiles, a distinct red asymmetry is observed at ra- 
dial velocities > 2000 km/s, the red wing being brighter 
than the blue one. This is indicative of stronger variability 
of the red wing as compared to the blue one during the 
monitoring period. 

On the rms Ha profile, a small peak is seen at a zero 
radial velocity, which is caused by an improper subtraction 
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Fig. 4. The H/3 (left) and Ha (right) profile decomposi- 
tion of the variable (top) and constant (bottom) compo- 
nents. The flux of the constant components is in units 
10~ 15 ergcm~ 2 s _1 A -1 . For the variable components, the 
relative increment of the line flux when the continuum 
flux increases by an amount of 1 x 10~ 15 ergcm~ 2 s _1 A -1 
is shown. The correlation coefficient between the variable 
component of lines and the continuum flux at different 
wavelengths along the line profile is shown in the middle 
panels. 




-1 0000 -5000 5000 10000 -1 0000 -5000 5000 1 0000 



Radial velocity (km/s) 

Fig. 5. The mean (top) and rms (bottom) H/3 (left) and 
Ha (right) profiles. The absissae are the radial velocities 
relative to the narrow components; the ordinates give the 
fluxes in units 10 -15 ergcm~ 2 s _1 A -1 . 



of the narrow Ha component for spectra of lower resolu- 
tion. The FWHM value of the mean and rms profiles is 
~ 6300 km/s and ~ 5800 km/s, respectively. These values 
are close to those obtained by Wandel et al. I|1999fl from 
spectra obtained during the 1989 - 1996. The mean pro- 
files of H/3 and Ha show blue peaks brighter than the red 
ones during the 1996 - 2001. 
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3.2.3. Time evolution of the broad H/3 and Ha 
profiles. 

The study of the shape of the broad emission line profiles 
and their variation in time, can help in choosing a suit- 
able model of the BLR. From our spectra it is seen that 
within every month, the broad line profile did not vary at 
a noticeable level. The flux in the lines also varied very 
slightly (as a rule, by 2 - 5%, and only in some cases up 
to 10%). Within every year, the shape of the profile did 
not vary at a noticeable level, but the fluxes in the broad 
emission lines varied considerably in 2000-2001. Then av- 
eraging the broad profiles over months or years, in order 
to increase the signal to noise ratios, allows us to see their 
time evolution in a more reliable way than on the individ- 
ual spectra. Therefore, for each month and each year, we 
have obtained the mean profiles. 

The averaged profiles of H/3 and Ha broad emission 
lines in subsequent years are shown in Fig. The evo- 
lution of the profile is well seen: in 1996, double peaks 
are present at radial velocities ~ —1000 km/s and ~ 
+ 1200 km/s; in 1997 - 1998, the double peaks are dis- 
tinctly seen at radial velocities of about ±1000 km/s, 
and in 1999, they are located at radial velocities of about 
±(500 - 900) km/s; in 2000 - 2002, a bright blue peak is 
still present at ~ -(700 - 1000) km/s, but on the red side 
at ~ +1000 km/s, a bending shoulder is seen instead of a 
bump. The relative brightness of the peaks varies: the blue 
peak was brighter than the red one in 1998 - 2002; but in 
1996, the red peak at a radial velocity of ~ +1200 km/s 
became brighter. These effects are best seen on the H/3 
profile than on the Ha one. It is because the narrow com- 
ponents are not well subtracted in the Ha case. 

In 2000 - 2002, a new bright bump is clearly seen in the 
red wing of the broad lines at a radial velocity of about 
~ +2500 km/s. In order to investigate whether or not 
the radial velocity of the new bump varies, we compared 
the observed profiles of H/3 obtained at different times 
(before subtraction of the narrow component). We have 
measured the radial velocity of the peak of this bump 
using good spectra with similar resolution (8-9) A and 
a high S/N ratio >50. We have also used the measure- 
ments of the peak location in May, 2003 obtained from 
the spectra taken with the 6 m telescope. The bump is 
relative broad, and the determination of the radial ve- 
locity of its peak is somewhat uncertain. Using the indi- 
vidual spectra of H/3 we defined the radial velocities of 
the bump peak by two methods: 1) as the mean-weighted 
location of the point in the bump corresponding to the 
level ~ 0.8Imax (Imax being the maximum of the peak 
intensity); 2) by fitting the bump top with a parabola. 
Both methods gave similar results. The only difference is 
that the second method (fitting by a parabola) gives ra- 
dial velocities of the bump peak systematically ~ 100 km 
higher, probably because of the slight blue asymmetry. 
We think that the first method gives more realistic results. 
Therefore we present the year-averaged radial velocities of 
this red peak derived by method 1 in Table El 1 — year; 2 
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Velocity(km/s) Velocity (km/s) 

Fig. 6. Annual averages from 1996 to 2002 of the observed 
profiles of H/3 broad emission component after subtraction 
of the continuum and the narrow components of H/3 and 
[Om]AA4959, 5007 (left). The same for Ha after subtrac- 
tion of the continuum and the narrow components of Ha, 
[N n] AA 6548, 6584 and [S n] AA 6717, 6731 (right). The ver- 
tical lines correspond to the radial velocities: dashed lines: 
±1000 km/s; dash-dot lines: ~ +2600 km/s; thin line: 
km/s. The absissae give the radial velocities relative to 
the narrow component of H/3 and Ha. The ordinates give 
the relative fluxes. The profiles are shifted vertically by a 
constant value. 



— the average radial velocity relative to the narrow com- 
ponent of H/3 and the root-mean-square error obtained 
from measurements on the individual spectra; 3 — the 
average continuum flux at the observed wavelength 5190 
A and its dispersion; 4 — the S/N ratio on the average 
spectrum in the region (5160-5220)1. One sees clearly on 
this table that the radial velocity of the peak decreases: in 
2000-2001 it corresponds to ~ (2500-2600) km/s, and in 
2002-2003 to ~ 2000 km/s, within the uncertainties. This 
effect is well seen on Fig. 7, where the year-averaged nor- 
malized profiles of H/3, derived by dividing the observed 
profiles by the average H/3 flux for each year, are presented. 
A similar result is obtained from the observed profiles af- 
ter subtraction of the narrow components. Thus, during 3 
years (from 2000-2001 to 2002-2003) a considerable shift 
~ 500 km/s of the red peak is observed on the H/3 profiles. 

From Table El there is no apparent correlation between 
the variations of the radial velocities of this peak and the 
average continuum flux. 
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Table 7. Year-averaged radial velocity of the blue and red bumps derived from the difference H/3 and Ha profiles 



Year 


H/3 


radial velocity (km/s) 


Ha radial velocity (km/s) 




v b 


Vrl 


v r2 


v bl 




v r2 


1996 


-906 ± 131 


+1461 ± 56 




-1200+ 69 


+1379 ± 103 




1997 


-833 ± 185 


+1216 ± 164 




- 696 + 32 


+1079+ 64 




1998 


-994 ± 91 


+1137+103 




-1075+ 67 


+1026+ 71 




1999 


-704 ± 200 


+ 838 ± 155 




- 944 ± 135 


+ 764 + 119 




2000 


-544 ± 134 


+817 : 


+2863 ± 262 


- 775 ± 192 


+800 : 


+2787 ± 261 


2001 


-704 ±316 


+954 : 


+2707+ 70 


-1043 ± 172 


+800 : 


+2526+ 92 


2002 


-938 ± 394 


+674 : 


+1898 ± 200 


-1158 + 144 


+584 : 


+1708 ± 191 



Table 6. Year- averaged radial velocity of the new red 
peak in 2000-2003, derived from the observed H/3 profiles 
before subtraction of the narrow component 



Year 


v r 


F(cont)* 


S/N 




(km/s) 


5190A 


(5160-5220)A 


2000 


2660+106 


9.27+1.28 


70 


2001 


2438+83 


9.66+1.26 


95 


2002 


2088+46 


7.11+0.98 


76 


2003 

± 7-1/ 


1986+75 


10.43 


95 

-2 I 'i -In 



* - F(cont) in units of 10" ib (ergcm"' i s" 1 A _i ) 



3.2.4. Ha and H/3 difference profiles 

For the analysis of the shape changes in the broad line 
profiles, the best method is that of direct subtraction of 
spectra, i.e. obtaining the difference spectra. In this case, 
the narrow emission lines and the absorption lines of the 
host galaxy, which may distort the observed emission line 
profile, are fully cancelled. In order to carry out a proper 
subtraction, it is necessary that the spectra had compara- 
ble spectral resolutions. To obtain difference spectra, one 
usually subtracts either the mean spectrum or a spec- 
trum representative of the minimum activity state. We 
have adopted the second alternative, as the spectrum in 
the minimum state, provides us with crucial information 
about the contribution of the narrow emission and ab- 
sorption lines, which are to be removed in order to carry 
a proper study of the profile shape changes with time. 
Hence, we obtained the H/3 and Ha difference profiles by 
subtraction of the minimum activity state spectrum, rep- 
resented by the mean spectrum of the May 15, 17 and June 
4, 2002 observations for H/3, and June 4, 2002 for Ha. The 
subtraction was performed by our scaling program (see 
Appendix A) , and a spectrum in minimum state was used 
as a reference. Since within every year, the shape of the 
line profiles changes only slightly, we have calculated the 
annual mean difference profiles for H/3 and Ha (Fig. (SJ. 
Variations of the profiles similar to those on the observed 
H/3 and Ha line profiles (see section 15. 2.3[1 are well seen. 
Peculiarities in the shape of the bumps in the difference 
profiles arc better seen than on the observed profiles. In 
the Ha difference profiles for 2000-2002, at radial velocity 
~ +1000 km/s a drop in brightness (dip) is seen, this is 




-10000 -5000 5000 10000 

Velocity [km s _1 ] 



Fig. 7. Year-averaged normalized profiles of H/3 in 2000- 
2003. The vertical dashed line corresponds to a radial ve- 
locity +2600 km/s. The absissae give the radial velocities 
relative to the narrow component of H/3. The profiles are 
shifted vertically by a constant value. 

due to over-subtraction of the [Nil] A6584 line. This due 
to the fact that the FWHM of this line is smaller than the 
one for the [OIII] A5007 line, used to reduce the spectra 
to similar resolution (see Section 2.4). This effect becomes 
noticeably on the 2000-2002 spectra, when the flux in the 
broad components decreases by factors of the order of 2 
to 3. From the monthly-averaged difference profiles of H/3 
and Ha, we defined the radial velocities of the peaks as 
the mean weighted location of the bump tops (for details 
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Fig. 8. Annual averages from 1996 to 2002 of the differ- 
ence profiles of the broad H/3 (left panel) and Ha (right 
panel) lines. The vertical lines correspond to the radial 
velocities: dashed lines: ±1000 km/s; dash-dotted lines: 
~ +3000 km/s; thin lines: km/s. The absissae represent 
the radial velocities relative to the narrow component of 
H/3 and Ha. The ordinates give the relative fluxes. The 
profiles are shifted vertically by a constant value. 

see Section 3.2.3). Our results are listed in Tabled there 
the annual averages for the radial velocities of the blue 
(Vf,), and red (V r i and V r 2) peaks along with their corre- 
sponding dispersion values are presented, for H/3 and Ha. 
In Table [7| it is seen that the peak radial velocities of H/3 
and Ha for different years agree with each other within the 
errors. It may be noted that during 1996-1999, the radial 
velocities of the blue (Vf,) and red bumps (V r i) were, on 
average, close to ±1000 km/s. However, it is distinctly 
seen that in 1996 the radial velocity of the red peak (V r i) 
was larger than in subsequent years. During 2000 - 2001, 
a decrease of the radial velocity of the blue (Vf,) peak in 
the H/3 difference profiles was observed. At that time the 
radial velocity determinations of the red peak (V^-i) from 
the difference profiles in 2001-2002 are uncertain (marked 
with a colon in Table [TJl, because in this velocity zone, 
a bending shoulder is present in the H/3 profile, and the 
previously mentioned dip in the Ha profile, due to over- 
subtraction of [Nil] A6584, affects our results. In 2000 a 
new distinct red bump appeared with a radial velocity of 
about +2800 km/s (V r2 in Table 0). By 2002 its radial 
velocity had decreased to about +1800 km/s, while its 
brightness became similar to that of the blue bump. We 
can definitely say that in the year 2000 a new bump in the 



red wing of the lines appeared and that its radial velocity 
decreased by about 1000 km/s between 2000 and 2002. 

3.3. Flux variability for different parts of Ha and Hf3 
profiles 

As mentioned above, the flux of the H/3 and Ha emis- 
sion lines as well as the AA5190 continuum flux, varied 
significantly between 1996 and 2002. Some parameters of 
the variability were noted in Section f3.ll However, as it is 
shown in Sections 13.2.31 13.2.41 important details such as 
bumps or oblique shoulders are present in the profiles of 
the broad lines. Therefore, it is of great interest to study 
the behaviour in time of these parts of the profiles relative 
to each other, and with respect to continuum variations. 
To this purpose, we have divided the emission line pro- 
files into several radial velocity bins and constructed the 
data sets for a number of time series that are of interest. 
We have chosen symmetrical bins relative to zero radial 
velocity of the H/3 and Ha line profiles. These bins in- 
clude either distinct peaks or notable features that were 
observed at different times (see Sections 13.2.31 l3~2.4[) : the 
H/31, Hal set from -3000 km/s to -2000 km/s; the H/32, 
Ha2 set from -1500 km/s to -500 km/s; the H/33, Ha3 
set from +500 km/s to +1500 km/s ; the H/34, Ha4 set 
from +2000 km/s to +3000 km/s. The flux associated 
to those radial velocity bins are named as: Fluxl, Flux2, 
Flux3, Flux4, respectively. Their light curves for H/3 and 
Ha are plotted in Figs. I9TTU1 ( left). 

The upper left and right panels present the light curves 
for the continuum near H/3 (Fig. 01 and Ha fFig. HU|) . From 
these figures, one can see that the flux in different parts 
of H/3 and Ha change in an almost identical manner. 

3.3.1. The shapes of the Ha and H/3 in relation to 
their variability 

The broad H/3 and Ha emission lines show changes in 
both: flux and shape of the profiles. Wanders & Peterson 
(1996) suggested a method for studying the variations of 
the profile shape by means of a normalization of the emis- 
sion lines by their total flux and then see how these nor- 
malized profiles change in the course of time. They named 
shape the function Fq(v,t), defined as: 

Fq(v, t) = F(v, t) - [< F(v, t) > / < F(t) >} ■ F(t), (5) 

where F(v,t) is the flux in the profile at a radial ve- 
locity v at a time t; < F(v, t) > and < F(t) > are the 
time-averaged (over the whole set of observations) flux at 
a radial velocity v and the total flux in the line, respec- 
tively; F(t) is the total flux in the emission line at a time 
t. By definition one gets J Fq(v,t)dv=0, if the integration 
is performed over the whole emission line. The shape value 
in a given radial velocity interval (or the shape value in 
a given profile part) is defined by integrating equation [5] 
in this velocity interval. From equation it follows that 
the value of Fq(v, t) will be close to zero when the profile 
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Fig. 9. Variations of the flux (left panel) and shapes 
(right panel) in different parts of the H/3 profile, and in 
the "combined" continuum at the observed wavelength 
A5190A. The radial velocity intervals of the profile seg- 
ments: the H/31, Hal set goes from —3000 km/s to 
-2000 km/s; the H/32, Ha2 set goes from -1500 km/s 
to —500 km/s; the H/33, Ha3 set goes from +500 km/s to 
+ 1500 km/s ; the H/34, Ha4 set goes from +2000 km/s 
to +3000 km/s. The flux in the lines are given in units 
of 10 _13 ergcm -2 s _1 , and in the continuum are given in 
units of 10~ 15 ergcm~ 2 s _1 A -1 . The shapes are in relative 
units (see Section l3.3.1fl . 
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Fig. 10. Same as Fig. El but for Ha and for the continuum 
at A6670. 



part F(v,t) changes proportionally to the total flux F(t). 
If the proportionality is violated due to delay effects or to 
some other reason, then Fq(v,t) will be either positive or 
negative. A positive value Fq(v,t) means that at a time 
t the given profile part is more prominent than in the av- 
erage profile, while a negative value means that this part 
is less prominent. So, F(v, i) is a sensitive indicator of the 
relative prominence of a part of the profile. If the profile 
shapes vary due to reverberation effects, Fq(v,t) must be 
correlated with the flux variations of the continuum after 
some delay. 

We have considered four shape time series in the ra- 
dial velocity intervals mentioned in section ^. 31 shapel for 
(-3000, -2000) km/s; shape2 for (-1500, -500) km/s; 
shape3 for (+500, +1500) km/s and shape4 for (+2000, 
+3000) km/s). 

The changes of the shape of the different profile parts 
are shown on the right panels of Figs. OH and for H/3 
and Ha, respectively. The changes in shape behave differ- 
ently from continuum variations: sometimes they change 
following the continuum changes, and some other time in 
quite the opposite manner. The behavior of the shape of 
different parts of the lines differs, though the changes of a 
given shape for Ha and H/3 are similar. 

3.3.2. The relations flux-shape and shape-continuum 
for the different parts of the H/3 and Ha broad 
emission lines 

A clear-cut distinction in the behavior of the flux and the 
shapes can be seen in Fig. 1111 where we show the flux- 
flux, shape-shape and shape-continuum flux correlations 
for the different parts of the emission profiles of H/3 and 
Ha. 

There is a high degree of correlation between the flux 
in the various parts of H/3 and Ha (Fig. El top)). The 
correlation coefficient values being 0.94 - 0.98. Yet, the 
correlations for the shapes of various parts of the lines are 
quite different. For example, 

— The shape4 and shape2 have a tendency to be anti- 
correlated with each other. 

— Shape4 in H/3 and Ha, including the bump at radial 
velocity +2500 km/s, does not correlate with the con- 
tinuum. 

— Shapc2 in H/3 is weakly related to the continuum, but 
the shape2 in Ha is not. 

— Shape3 and shape2 have a tendency to be correlated 
in H/3 and weakly anti-correlated in Ha. 

— A good correlation with the continuum for the shapc3 
in H/3 and Ha is observed. 

Note that the analysis of shape2 and shape3 in Ha, 
should be done with care, because the narrow components 
of Ha and nitrogen are not properly subtracted in some 
profiles. Hence, the derived fluxes are subject to uncertain- 
ties, and the profiles in these regions may be considerably 
distorted. In the H/3 case, the narrow component is well 
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Fig. 11. Flux-flux, shape-shape and continuum flux - 
shape correlations for different parts of the Ha and the 
H/3 lines (the number shape is inside the plot). The cor- 
relation coefficients r are given inside the plot. 

subtracted and data are more reliable. Therefore, from the 
H/3 analysis alone, for shape2 and shape3 that include the 
double-peaks at radial velocities close to ±1000 km/s, we 
conclude that they are weakly correlated with each other 
and also with the continuum. 

Thus, the shape of the individual parts of the line pro- 
files changes, unlike their fluxes, in a more complicated 
way. So, since the changes in shapes correlate sometimes 
slightly with the continuum changes, or do not correla- 
tion with them, then the changes may not be related to 
reverberation effects, and some other mechanisms should 
be invoked for their explanation. 

3.4. The Balmer Decrement 
3.4.1. Integral Balmer Decrement 

^From the spectra with the subtracted narrow component 
template (see section |2~4"|) . we determined the fluxes of Ha 
and H/3 emission line broad components within the radial 
velocity intervals ±6000 km/s. The integral Balmer decre- 
ment (BD) is the flux ratio F(Hct)/F(H0). Note that the 
Fell (42 multiplet) line at the rest wavelength 4924 A cor- 
responds to a radial velocity of +3874 km/s relative to 
the narrow component of H/3. This line sets in the broad 
red wing of H/3 close to the blue wing of the emission 
line [OIII] 4959. No measurable feature is seen on our 
spectra at this place. Since it is known that the inten- 
sities of FeII42 (4924, 5018, 5169) lines are comparable, 
we have estimated the contribution of Fell42 5018 A line. 
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Fig. 12. Variations of the Balmer decrement 
F(Hot)/F(H(3) (top) and of the continuum flux (middle) 
in 1996 - 2001; variations of the Balmer decrement versus 
those of the continuum flux (bottom). 



This line is well seen on our spectra. We have obtained 
for the ratio of FeII42(5018) to the flux integrated within 
the interval of radial velocities ±-6000 km/s (i.e. onto the 
broad component of H/3): F(FeII42)/F(H/3) - (0.02-0.03) 
in 1996-1999, and - (0.03 - 0.06) in 2000-2001. Thus, 
the contribution of Fell42 to the H/3 broad component is 
- (2 - 3)% in 1996-1999, and - (3 - 6)% in 2000-2001. 
This is within the measurement errors of the Balmer decre- 
ment. Therefore, we did not take into account the contri- 
bution of FelI42. Fig. El shows the behavior of Balmer 
decrement (upper panel), the A 5190 continuum, (middle 
panel), and the relation between changes of the Balmer 
decrement and continuum flux (bottom panel). There an 
anticorrelation between the changes of the Balmer decre- 
ment and the continuum flux is evident (correlation coef- 
ficient ~ 0.52, bottom panel in Fig. I12J) . 

The same tendency is present in the annual averaged 
values of the BD and the continuum fluxes. In Tabled the 
mean values of the Balmer decrement are listed for several 
years together with the mean values of the continuum flux. 



Shapovalova, A.I., et al.: H/3 Ha monitoring of NGC 5548 



13 




1996 Feb. 13 




2 - 2000 Feb, 26 




2 - 2001 Mar, 12 



-4000 -2000 



2001 Jul, 12 



2000 4001 -4000 -2000 




1998 Feb, 22 




During 1996 - 1999 in the radial velocity interval from 
~ -4000 km/s to ~ +2000 km/s, the Balmer decrement 
varied only slightly and had a value of ~ 4.0, while in the 
radial velocity interval +(2000 — 4000) km/s, it decreased, 
from ~ 4.0 to 2.0. In the year 2000, the Balmer decrement 
increased, as a whole. Yet, a monotonic decrease with ra- 
dial velocity is clearly present. During 2001, maximum 
values (<~ 5) were observed at the line center, and at the 
edges the BD decreases with a somewhat larger gradient 
in the red wing. 

4. Summary of our main results and comparison 
with other results from the literature 



- 2000 Apr, 25 




2000 400C 



Fig. 13. Variations of the Balmer decrement 
F (Ha) / F (H.p) along some individual profiles of the 
broad lines. The absissae represents the radial velocities 
relative to the narrow components. 

Table 8. The anuual averaged Balmer decrement and con- 
tinuum fluxes F(cont) at the observed wavelength A5190 A 



Year F(Ua)/F(U/3) 



F(conty 



1996 
1998 
1999 
2000 
2001 



3.48+ 0.45 
3.51+ 0.44 
3.51+ 0.44 
4.22± 0.42 
4.33+ 0.58 



F(cont) in units of 10" 15 (erg cm s 



14.62+ 2.2 
16.96+ 1.69 
14.05+ 1.62 
9.27+ 1.28 
9.66+ 1.26 



iFvom Table H it is seen that during 1996 - 1999 the 
average flux in the continuum varied slightly, while the 
Balmer decrement did not changed. However, in 2000 - 
2001 the average flux in the continuum decreased and the 
Balmer decrement became considerably larger (steeper). 



3.4.2. Balmer Decrement across some individual 
profiles 

Fig. 1131 shows the Balmer decrement across some individ- 
ual broad lines profiles. The decrement is similar to that of 
other profiles for the same years. Yet, the Balmer decre- 
ment across the line profile changes from year to year. 



1. 



In Section E3 we have shown that the observed mean 
and rms profiles as well as the annual averaged, ob- 
served and difference profiles of Ha and H/3 present a 
double-peak structure, i.e. bumps at radial velocities 
~ ±1000 km/s ( zero velocity being defined as that of 
the narrow components) . While the double peaks were 
obtained after subtraction of the narrow line compo- 
nents, a question arises about the reality of these fea- 
tures, as it dependent on the method of subtraction 
of those components. The following facts are taken as 
evidence to the reality of the double peaks: 1) Double- 
peaks are obtained on the variable component when 
the original H/3 and Ha profiles are decomposed, be- 
fore subtraction of the narrow components (Section 



2) The appearance of double peaks on the dif- 
ference profiles of H/3 and Ha, when the narrow lines 
are subtracted automatically fSection |3.2.4|l . 
The relative intensity of these peaks varied: in 1996, 
the red peak was brighter than the blue one, and vicev- 
ersa, in 1998-2002, the blue peak became brighter. We 
have seen that the radial velocities of the double peak 
vary within - ±(500 - 1200) km/s. 
The double-peaked structure in the central part of the 
broad lines in the radial velocity range ~ ± 500 - 
1500 km/s was observed earlier. Double-peaks at ve- 
locities ~ ±500 km/s were seen in the annual aver- 
aged profiles of H/3 in 1986 - 1987 by Wanders and 
Peterson (|1996j) . On the individual difference profiles of 
H/3, the double peaks were also evident in 1992 (Iijima 
& Rafanelli 1993). Very distinct double peaks were also 
noticeable on the difference Ha profiles in July 1986 - 
June 1985 (Stirpe et al. rOM| and in May 1987 - July 
1986 (Stirpe & de Bruyn HMD). 

During 2000 - 2002, when the flux of the continuum 
and the lines greatly decreased, a new distinct peak 
at a radial velocity of ~ +2500 km/s appeared on 
the profiles of Ha and H/3. This peak had compara- 
ble brightness to the blue peak at a radial velocity 
of < -1000 km/s)(Section A radial velocity 

decrease of ~ 500 km/s for this feature occured be- 
tween 2000 and 2002 (Table Eland Fig. 0). The same 
peak is seen on the difference profiles of Ha and H/3 
(Section I3.2.4J) . Its radial velocity decreased as well, 
but by a larger value (~ 1000 km/s) than the one de- 
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rived from the observed profiles. This is a result of the 
details present in the line profiles at minimum activity 
state, which are subtracted from the observed profiles. 
In the red wing of H/3 at minimum activity state, a 
distinct wide asymmetric weak bump is observed at a 
radial velocity of about +2200 km/s. Therefore on the 
difference profiles in the red wing, we observe excess 
emission relative to this weak bump in the spectrum 
at minimum state. Thus, from our data set, we can es- 
tablish the appearance of a new peak in the red wings 
of Ha and H/3 during 2000 - 2002. The radial velocity 
of this peak in 2002 decreased by ~(500 - 1000) km/s, 
with respect to the value determined for the year 2000. 

3. /From the decomposition of Ha and H/3 profiles (see 
Section 13.2. 1|) we have found that the fluxes of the 
variable component of the lines correlate well with the 
flux variations of the continuum. In the radial velocity 
interval ~ —4000 km/s to ~ +5000 km/s the corre- 
lation coefficient has a high value r ~(0.8 - 0.9). In 
Section 13.31 we have shown that the flux of different 
parts of the profiles, change in an almost identical man- 
ner being highly correlated both with each other (r ~ 
0.94 - 0.98) and with the continuum flux levels (r ~ 
0.88 - 0.97). These facts indicate that the flux variabil- 
ity in different parts of the line profiles on short time 
scales is caused mainly by the reverberation effect. 

4. Line profile shapes evolve slowly, on time scales of 
months to years. The shape variations of the different 
parts of the line profiles, mildly correlate with each 
other and/or with the continuum variations, or simply 
do not correlate at all (see Sec. 13.3.21 and Fig. Ill|) . So, 
the profile shape changes on long time scales are not 
due to reverberation, as earlier discussed by Wanders 
& Peterson (|T§9l)|> . 

5. The flux ratio of F(Ha)/F(H/3) (BD) shows a clear 
anti-correlation with the continuum flux (see Fig. El 
and Table |HJ. Such an effect has also been observed in 
other monitoring campaigns of Sy galaxies (Dietrich et 
al. 119931 Shuder^981). In addition, there are changes 
of the Balmer decrement across the profiles and it 
varies from year to year (see Fig. 113(1 . 

5. Discussion 

The study of the Ha and H/3 broad line profiles in NGC 
5548 with different techniques, shows the presence of 
double-peak structures at radial velocities ~ ±1000 km/s 
relative to the narrow steady component (Section point 

I)- ^ 

The appearance of double peaked structures in the pro- 
files of broad emission lines is predicted by some types 
of models: a) different versions of accretion-disc models 
(Dumont & Collin-Souffrin I199()al I1990bl Rokaki et al. 
113321 Chen & Halpern 1989, Eracleous et al. 11335) ; b) bi- 
conic gas flows (Veilleux and Zneng lT991[l : c) binary black 
holes (Gaskell 1X353) . 

The strongest argument against the binary black hole 
model is the fact that the velocity curve of the Ha red 



peak in 3C390.3 (the most prominent object with double 
peaks) is inconsistent with the best-fitting binary mod- 
els (Eracleous et al., 1997). Within the framework of the 
model of binary black hole, the appearance of the third 
peak in NGC 5548 at a radial velocity of ~ +2500 km/s 
is impossible to explain. The biconical gas flows imply 
mainly radial motions in BLR, not supported by any ob- 
servations. For instance, from the CCF analysis of the H/3 
emission line profiles from the data of AGN- Watch mon- 
itoring in 1989-1993, Wanders and Peterson <|1996|l ex- 
cluded the possibility of large radial motions in BLR of 
NGC 5548. Kollatschny and Dietrich (JT996) showed that 
the outer wings of the Balmer lines have a tendency to re- 
pond faster than their cores to the continuum variations. 
We also tried to carry out the CCF analysis for different 
parts of the Ha and H/3 profiles using our data. But be- 
cause of our poor sample there were large uncertainties in 
the time delay determination. However, in spite of this, we 
noted that Ha and H/3 show a same tendency for smaller 
time delays at high velocities delays than at low veloci- 
ties, and that there is no time delay between the blue and 
red wings of H/3 within 3-4 days. Of course, our results 
are only qualitative, but they do not contradict the re- 
sults of the AGN- Watch monitoring. Thus, the available 
observations indicate predominantly rotational (or virial) 
motions of NGC 5548 in BLR, and support the scenario of 
formation of the broad emission lines in an accretion disk. 
Assuming a standard geometrically thin disc, Dumont and 
Collin-Souffrin ( 1990ai n"990b|l have shown that the broad 
Balmer lines can be almost entirely due to disc emission. 
They obtained a large variety of profiles, including dou- 
ble peak symmetrical profiles for non-relativistic cases. 
Relativistic corrections introduce profile asymmetries: the 
blue side of the profile becomes brighter than the red one. 
However, in 1996, the red peak was brighter than the blue 
one, while at other times (1998-2002), the blue peak was 
the brighter one. The cases when the red peak is brighter 
than the blue one, are in contradiction with the predic- 
tions of relativistic circular disc models. However, as the 
calculations of Eracleous et al. Q1995[l show, such case is 
possible, when the broad lines originate in a relativistic, 
eccentric disc. In this case the changes in the line profile 
shapes, will be mainly due to precession of the disc. If this 
is the case, the radial velocities and the intensities of the 
double peaks would change slowly (in timescales of months 
to years), without following the effects of reverberation. 

An alternative explanation for the observed intensity 
variations of the double peaks relative to each other, is 
to involve local disc inhomogeneities, which may be re- 
sponsible for some of the substructures of the emission 
line profiles. Actually it is now clear that at the BLR dis- 
tance, the disc becomes gravitationally unstable (Collin & 
Hure 1993J). At that point the disc is expected to be highly 
inhomogeneous, including clouds and possibly spiral arms 
and it is certainly not constituted by a uniform density 
medium (e.g. Bunk, Livio & Verbunt 1 991) Chakrabarti 
& WiitaUHMI). 
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Another issue which must be taken into consideration, 
is the fact that the line emitting regions should be illumi- 
nated by the central continuum source, in order to be able 
to reprocess an important fraction of such radiation. It is 
now well established that the UV-X source has small di- 
mensions as compared to the BLR (say 10-100i?Gi where 
Rc is the gravitational radius, the BLR being located at 
10 3 ~ 4 i?c)- So, if the BLR constitutes itself the outskirts 
of the disc, it should be able to see the central source. 
There to this aim, there are several possibilities: 

— As suggested by Dumont & Collin (1990a, 1990b), the 
central radiation can be back-scattered towards the 
disc, as most of the hard radiation is emitted in the soft 
X-ray range (<s 10 keV). Yet, this requires a highly ion- 
ized Thomson thick medium, surrounding the disc and 
the central source. This medium would create a low en- 
ergy cut-off (at many keV), which is not observed in 
the high energy spectrum of AGNs and would change 
dramatically the hard X-range spectrum due to multi- 
ple Compton scattering. 

— The disc should extend vertically, up to a scale height 
equal to at least one tenth of the radius. This would im- 
ply that a good fraction of the observed velocity (tens 
of percent) shall be non-rotational in origin. Probably 
due to large scale chaotic gas dynamical motions - 
macro-turbulence- added on top of regular rotation. 

— The disc can be geometrically thin but warped, as sug- 
gested by recent observations (Kinney et al. 2000). The 
warping could be due to self-irradiation instabilities 
(Pringle 1996), to misalignment of the spin of the BH 
with the disc (Bardeen & Petterson 1975), or to disc 
feeding through a misaligned inflow. 

In 2000 - 2002, we detected a new peak at a ra- 
dial velocity of ~ +2500 km/s in the red wing of Ha 
and H/3. This bump moved gradually along the line pro- 
file and its radial velocity decreased between 2000 and 
2002 by ~ 500 km/s on the observed profiles and by 
~ 1000 km/s on the differential profiles (Section QJ. 
On the normalized H/3 profiles (1989 - 1993), Wanders 
and Peterson (119961 noted the periodical appearance of 
broad red- and blue- wing shoulders at radial velocities of 
~ ±2500 km/s. According to these authors, their relative 
intensities change over the years, but not their positions in 
the line profile. They interpreted this fact with a model of 
cloud illumination by an anisotropic source. However, on 
the normalized H/3 profiles shown by these authors, no dis- 
tinct peaks (bumps) are noted. It looks more as an asym- 
metry in the line profiles, i.e. a brightness increase within 
the velocity interval ±1000 — 2500 km/s. It was impossi- 
ble to determine the location of the peaks from these pro- 
files, since they were absent. Possibly, during the period of 
1989-1993, the components appearing periodically in the 
velocity interval ±(1000—2500) km/s were very broad and 
of low contrast. Therefore, the averaged profile obtained 
after summation of a large amount of spectra, showed the 
brightness increase in the form of inclined shoulders. Note 



that in 1991-1993 Sergeev et al. (1994) observed a dis- 
tinct low contrast broad emission feature in the far red 
wing of the Ha profile. The radial velocity of this feature 
decreased from ~ +4500 km/s in 1991 to ~ +2100 km/s 
in 1993. So, from our observations and those of Sergeev 
et al. I|1994p . we deduce that the position of the bump 
changes with time in the line profile, and therefore the 
interpretation of Wanders and Peterson {H96) based on 
a constant radial velocity of the bump is not correct. We 
think that the appearance of such details (bumps) moving 
along the profile is due to inhomogeneities in the disk. For 
time scales of years, which correspond to the dynamical 
time of the BLR, one can invoke to account for these in- 
homogeneities, a structure like a spiral wave at ~ 10 light 
days from the center which would rotate at approximately 
the same speed as the disc, causing a shift of the peak in 
the profile in from 1000 to 2500 km/s in a few years (for 
a black hole mass of 7 x 10 7 Mo, the rotation timescale is 
on the order of 20 years) . Such a single spiral arm is also 
invoked to explain the profiles and their variations in the 
LINER/Seyfert 1 galaxy NGC 1097 (Storchi-Bergman et 
al. 2003). A spiral wave will make a "bump" in the disc, 
as it thicker than the rest of the disc. Therefore it would 
be more illuminated than the rest of the disc, causing a 
localized enhancement in the line profile. The rotation of 
the spiral wave should be accompanied by a change of the 
time lag between the continuum and the lines. In this case 
there is no obvious link between a decrease of the contin- 
uum flux level and a velocity shift of the bump. 

There are however, other possible explanations differ- 
ent from that of a spiral arm. For instance one can in- 
voke a hot spot due to a collision of the disc with a pass- 
ing star (e.g. Syer, Clarke & Rees 119911 Chakrabarti & 
Wiita rK)93f) . or the explosion of a supernova releasing a 
rcmnant(a "SN cloud") (Collin & Zahn 1999). In this case, 
the emitting material should be receding, with a velocity 
increasing with time. In the model of the SN cloud, one 
can explain the correlated decrease of the continuum by 
assuming that as the SN cloud grows, it begins to cover 
the line of sight of the central source (in such a dense 
environment, a supernova would be denser, with a much 
larger column density than a normal supernova, during a 
few years). In both cases the new material will be incor- 
porated later on to the inner disc and the accretion rate 
will increase. But this would take place a long time after 
the explosion of the supernova or the capture of the star 
by the disc, in the viscous timescale of the inner disc, i.e. 
in tens of years, so we cannot check this hypothesis at 
present. 

For the two peak structure, one could invoke a double- 
armed spiral shock structure as proposed by Chakrabarti 
and Wiita (fT5Mjl . 

In a detailed physical study of the Balmer line emitting 
region in NGC 5548, Dumont et al. (1998) considered a 
sample of observations spanning a relatively small time (9 
months) in 1989. They concluded that it must be heated 
partially by a non radiative mechanism, in particular be- 
cause the variations of the lines were too small compared 
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to those of the UV continuum in short term scales. Indeed, 
we noticed in the present paper that the line flux varia- 
tions are very small in short time scales. This non radia- 
tive heating can be due for instance to sound waves. Such 
waves can be easily provided by a spiral wave or by a 
supernova, or by a star colliding with the disc. 

So, we can suggest that the BLR is a highly turbulent 
medium in rotation constituting a kind of semi-thick inho- 
mogeneous disc. It is partly heated by the central UV-X 
source, and partly mechanically through the dissipation 
of turbulent motions. At short time scales, instabilities of 
the inner disc create rapid changes of the continuum flux, 
inducing the reverberation effects, but smoothed by the 
non radiative part of the heating which varies much more 
slowly. At longer time scales, structural changes of the 
BLR induce changes of the line profiles and of the rever- 
berated flux, and should lead after some years to changes 
in the accretion rate and consequently changes in the con- 
tinuum flux. The inhomogcneitics in the disc will boost 
the emissivity at specific radial velocities, thus producing 
bumps or asymmetries in the line profiles. However, one 
cannot choose among a variety of models of inhomoge- 
neous disc the one which suits better the observed evolu- 
tion of the Balmer line profiles in NGC 5548 in 1996-2002. 

The variations of the integral Balmer decrement can 
also be interpreted in this framework. In Section 0] (point 
5), we pointed out that the integrated Balmer decrement 
is anticorrelated with the continuum flux level. As the line 
flux in Ha and H/3 correlates well with that of the contin- 
uum, we infer that the change of the integrated Balmer 
decrement is also caused by changes in the continuum flux. 
Indeed, when the ionization parameter decreases for a con- 
stant density plasma, an increase of F(Ha)/F(H/3) inten- 
sity ratio is expected, (cf. for instance Wills et al. 1985). 
This is due to the decrease of the excitation state of the 
ionized gas: the temperature of the ionized zone being 
smaller, the population of the upper levels with respect 
to the lower ones decreases. However the great differences 
observed in the behaviour of the Balmer decrement along 
the individual profiles are certainly not due to the vari- 
ations of the ionizing flux, as they take place in longer 
time scales. They are more probably due to local density 
fluctuations of the inhomogeneous medium, leading also 
to variations of the ionization parameter, and inducing 
changes of the Balmer decrement across the line profiles. 

6. Conclusions 

Between 1996 and 2002, we have carried out a spectral 
monitoring of the Seyfert galaxy NGC 5548. It is found 
that the flux in the lines and the continuum gradually de- 
creased and reached minimum values in May- June 2002, 
i.e. the continuum source gradually faded. The line wings 
at maximum light states (1998) correspond to a Syl type, 
while at minimum light states (2002), they are similar 
to a Syl. 8 type. It was shown that the observed mean 
and rms line profiles, the variable part of the H/3 and Ha 
profiles, and the difference profiles present double peaks 



at radial velocities ~ ±1000 km/s. During 1996, the red 
peak was brighter than the blue one, and in other years 
(1998 - 2002) the blue peak became brighter than the red 
one. In 2000 - 2002, we observed a third distinct peak 
in the red wings of Ha and H/3, at a radial velocity of 
~ +2500 km/s. The radial velocity of this peak decreased 
from 2000 to 2002 by ~ 500 km/s on the observed pro- 
files, i.e. this peak moved gradually across the line profile 
towards lower radial velocities. The observed fluxes in dif- 
ferent parts of the broad lines vary quasi-simultaneously. 
The flux in different parts of the profiles, are well corre- 
lated with each other, and with the continuum flux as well. 
Yet, the change in shapes of the variuos parts of the lines, 
either they weakly correlate, or simply do not correlate at 
all with the changes in the continuum flux. Hence, these 
changes are not caused by reverberation effects. We have 
observed a general increase of the integral Balmer decre- 
ment with a continuum flux decrease. This probably due 
to a decrement in the ionization parameter, caused by a 
decrease of the flux in the continuum. The behavior of the 
Balmer decrement across the line profiles, differs greatly 
from time to time, is probably due to variations of density 
in inhomogeneities in an accretion disc. 

Section [5] tries to integrate all these observations in 
the framework of a model, and presents arguments in favor 
of the formation of the broad Balmer lines in a turbulent, 
partly mechanically heated accretion disc including large, 
moving "thick" inhomogeneities, capable of reprocessing 
the central source continuum. 
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Appendix A: Modification of the spectrum scaling 
method 

In the monitoring programs of AGN, the spectral scaling 
scheme described by van Groningen & Wanders ( 1992J) is 
used. The main idea of the algorithm is the creation of 
the difference spectrum between an input spectrum and 
a reference spectrum, for which the flux is assumed to be 
constant. The difference spectrum is represented by the 
simple analytical function (usually by a 2nd order poly- 
nomial). Then a \ 2 of this correspondence is minimized 
by a grid search method by successive variations of 3 in- 
put parameters: a flux scaling factor, a wavelength shift, 
and a difference in resolution of the spectrum (A FWHM). 
For the latter, a convolution with Gaussian function is se- 
lected. For data with the S/N ratio ~ 10, the error of 
the flux scaling factor is < 5%. However, the method 
of van Groningen & Wanders Q1992[l is unstable to the 
selection of initial parameters ( the zero approximation 
is done manually). In order to circumvent this problem, 
we have modified the method. The difference between the 
individual spectra (obj) and the "reference" one (ref) is 
represented by a 3rd degree polynomial, and for the mini- 
mization of the differences, a downhill simplex method by 
Nelder and Mead fl!965fl is used. The latter is more stable 
than the grid search method used by van Groningen & 
Wanders <|1992fl . As a zero approximation, the flux in the 
spectrum lines was determined automatically after sub- 
traction of a linear continuum determined by the begin- 
ning and the end of a given spectral interval. The scaling 
procedure is then carried out with the program means in 
IDL, the program is fast and stable. The program output 
is similar to that of van Groningen & Wanders (1992): the 
flux scaling factor, relative wavelength shift and Gaussian 
width which is used for convolution with one of the spec- 
tra for spectral resolution correction, values of x 2 an( i a 
for the power approximation to the spectrum difference, 
the scaled and difference spectra (obj-ref). The latter are 
obtained after reducing the spectra to the same spectral 
resolution. In order to check the correctness of the scaling 
method several tests have been carried out. 

1. We have tested the accuracy in the determination of 
the flux scaling factor by means of a model spectrum 
of different S/N ratios. As model data, we adopted 
the result of a multi Gaussian approximation to the 
Jan 21,1998 NGC 5548 spectrum, for the wavelength 
interval (4700 - 5400) AA including the [O in] AA 4959, 
5007 and H/3 emission lines. 

In testing, the oxygen line intensities remained con- 
stant while the H/3 intensity was varied between 10 
and 500 percent (H/3 = 100% corresponded to the ob- 
served spectrum of NGC 5548 in Jan 21 1998). A con- 
tinuum level was added to the lines and the spectrum 
was convolved with grey noise, to a chosen value of the 
S/N ratio. 200 experiments were carried out for each 
of the selected S/N ratios. From the simulated spectra 
with S/N=20, the average values of the flux scaling fac- 
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tor show systematic decrease from 1.015 to 0.995 for 
changes of the H/3 intensity from 10% to 500%. The 
dispersion of the flux scaling factor being about 2.5% 
As the S/N ratio is increased to 40 (typical value for 
spectra obtained in our monitoring campaign), system- 
atic errors and dispersion values went down to 0.5% 
and 1% respectively. 

2. With the model spectra we tested the method accuracy 
for the determination of the flux scaling factor in the 
case in which spectra of different resolution are com- 
pared. In this case the model spectrum adopted is the 
same as in case 1., with H/3 intensity 100%, but before 
introducing noise, a convolution with a Gaussian of the 
proper width is done. The value of the Gaussian disper- 
sion (ct=FWHM/2.35) varied from 0.01 to 5.59 pixels 
with steps of 0.01 pixel (dispersion being 2 A/px). We 
found that for spectra with S/N ratio ~ 20, the value 
of (7 > 0.51 and the flux scaling factor varies by less 
than 1%. We also checked for changes of the flux scal- 
ing factor in the observed blue spectra of NGC 5548. 
For this purpose, we first scaled the spectra adopting a 
spectrum with a resolution (~ 8A) as a referenced one. 
Then, the scaled spectra were processed by the pro- 
gram adopting a spectrum of lower resolution (~ 15A) 
as a the reference one. The values for the flux scal- 
ing factors so derived, coincided within an accuracy of 
~ 1%. 

3. The wavelength shift (AA) is always recovered with 
high accuracy by the method. 

Thus by scaling our AGN spectra by means of the 
modified method of van Groningena & Wanders (1992) 
as described above, one can obtain correct values of the 
scale parameters, their errors being dependent only on 
the quality of the spectra. 

Appendix B: The tables presented in electronic 
form 
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Table B.2. Log of the spectroscopic observations. Columns: 1 - number; 2 - UT date; 3 - Julian date; 4- code according 
to Table 1; 5 - projected spectrograph entrance apertures; 6 - wavelength range covered; 7 - spectral resolution; 8 - 
mean seeing; 9 - position angle (PA) in degrees; 10 - signal to noise ratio in the continuum (5160 - 5220) AA near H/3 
and (6940 - 7040) AA near Ha. 



No 


UT-Date 


JD 


Code 


Aperture 


Sp. range 


Res. 


Seeing 


PA 


S/N 






(2400000+) 




(arcsec) 


(AA) 


(A) 


arcsec 


(deg) 


5160-5220 AA 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


1996Janl4 


50097.570 


LI 


2.4x13.5 


3600-7200 


9 


4 




28 


2 


1996Janl5 


50098.581 


LI 


4.2x13.5 


3600-7200 


9 


4 




17 


3 


1996Febl4 


50127.579 


L 


1.5x6.0 


3100-5800 


6 


2 




95 


4 


1996Febl4 


50127.569 


L 


1.5x6.0 


4500-7200 


6 


2 




69 


5 


1996Febl4 


50128.432 


L 


1.5x6.0 


3100-5700 


6 


3 




114 


6 


1996Febl4 


50128.472 


L 


1.5x6.0 


4500-7200 


6 


2 




138 


7 


1996Marl9 


50162.382 


L 


2.0x6.0 


3600-5600 


8 


3 


40 


83 


8 


1996Marl9 


50162.399 


L 


2.0x6.0 


4700-7400 


6 


3 


40 


92 


9 


1996Mar21 


50164.395 


L 


2.0x6.0 


3600-5600 


8 


3 


45 


79 


10 


1996Mar21 


50164.412 


L 


2.0x6.0 


4700-7400 


6 


3 


45 


77 


11 


1996JullO 


50275.283 


L 


2.0x6.0 


3600-5400 


8 


1.4 


82 


96 


12 


1996JullO 


50275.353 


L 


2.0x6.0 


5200-7000 


8 


1.4 


82 


114 


13 


1997Apr05 


50544.414 


LI 


4.2x19.8 


4100-5900 


9 


4 


90 


40 


14 


1997Apr05 


50544.455 


LI 


4.2x19.8 


5500-7300 


9 


4 


90 


29 


15 


1997Apr08 


50547.331 


L 


2.0x6.0 


4444-5244 


4.5 


3 


6 


30 


16 


1997Aprl4 


50553.322 


L 


2.0x6.0 


4400-5300 


4.5 


4 


9 


41 


17 


1997Aprl4 


50553.330 


L 


2.0x6.0 


6200-7100 


4.5 


4 


9 


44 


18 


1998Jan21 


50834.632 


L 


2.0x6.0 


3800-6200 


8 


3 


62 


97 


19 


1998Jan21 


50834.639 


L 


2.0x6.0 


5900-7598 


8 


3 


69 


126 


20 


1998Feb23 


50867.535 


L 


2.0x6.0 


3800-6200 


8 


2 


53 


73 


21 


1998Feb23 


50867.539 


L 


2.0x6.0 


5900-7598 


8 


2 


58 


94 


22 


1998Apr27 


50931.474 


LI 


8.0x19.8 


5500-7300 


9 


4 




28 


23 


1998Apr30 


50934.438 


LI 


8.0x19.8 


5500-7300 


9 


4 


146 


48 


24 


1998Apr30 


50934.467 


LI 


8.0x19.8 


4000-5800 


9 


4 


146 


73 


25 


1998May04 


50938.293 


L 


2.0x6.0 


3700-6200 


8 


2.5 


22 


63 


26 


1998May04 


50938.289 


L 


2.0x6.0 


5900-7598 


8 


2.5 


20 


49 


27 


1998May06 


50940.433 


L 


2.0x6.0 


3700-6200 


8 


3 


108 


62 


28 


1998May06 


50940.436 


L 


2.0x6.0 


5900-7598 


8 


3 


107 


54 


29 


1998May07 


50941.551 


L 


2.0x6.0 


3700-6200 


8 


3 


93 


106 


30 


1998May08 


50942.452 


L 


2.0x6.0 


4100-5450 


4.5 


2 


104 


41 


31 


1998May08 


50942.461 


L 


2.0x6.0 


3700-6200 


8 


2 


102 


78 


32 


1998May08 


50942.472 


L 


2.0x6.0 


5700-7700 


8 


2 


100 


22 


33 


1998Junl7 


50982.407 


LI 


8.0x19.8 


4050-5850 


9 


2 





74 


34 


1998Junl9 


50984.287 


LI 


8.0x19.8 


4050-6147 


9 


4 





81 


35 


1998Junl9 


50984.328 


LI 


8.0x19.8 


5500-7300 


9 


4 





47 


36 


1998Jun20 


50985.350 


LI 


8.0x19.8 


4050-6000 


9 


2 





67 


37 


1998Jun20 


50985.281 


LI 


8.0x19.8 


5500-7300 


9 


2 





32 


38 


1998Jun25 


50990.409 


L 


2.0x6.0 


3600-6100 


8 


3 


84 


61 


39 


1998Jun26 


50991.331 


L 


2.0x6.0 


3600-6100 


8 


2 


102 


40 


40 


1998Jull4 


51008.778 


GH 


2.5x6.0 


3979-7230 


15 


2 


90 


106 


41 


1998Jull7 


51011.683 


GH 


2.5x6.0 


4210-7471 


15 


2 


90 


84 


42 


1998Jul26 


51020.703 


GH 


2.5x6.0 


3960-7231 


15 


2.2 


90 


73 


43 


1998Jul27 


51021.712 


GH 


2.5x6.0 


3930-7219 


15 


2.5 


90 


85 
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No 


UT-Datc 


JD 


Code 


Aperture 


Sp. range 


Res. 


Seeing 


PA 


S/N 






(2400000+) 




(arcsec) 


(AA) 


(A) 


arcsec 


(deg) 


5160-5220 AA 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


44 


1999Janl3 


51191.952 


GH 


2.5x6.0 


4150-7436 


15 


2.4 


90 


94 


45 


1999Janl4 


51193.011 


GH 


2.5x6.0 


4151-7438 


15 


2 


90 


94 


46 


1999Jan22 


51200.542 


LI 


4.2x19.8 


4000-5800 


9 


3 


90 


62 


47 


1999Jan23 


51201.579 


LI 


4.2x19.8 


5550-7350 


9 


3 


90 


57 


48 


1999Jan25 


51203.565 


LI 


4.2x19.8 


4000-5800 


9 


4 


90 


49 


49 


1999Jan26 


51204.576 


LI 


4.2x19.8 


5550-7350 


9 


4 


90 


59 


50 


1999Marl5 


51252.970 


GH 


2.5x6.0 


4200-7500 


15 


2 


90 


138 


51 


1999Mar24 


51262.444 


LI 


4.2x19.8 


4050-5850 


9 


3 


90 


68 


52 


1999Mar24 


51262.452 


LI 


4.2x19.8 


5550-7350 


9 


3 


90 


15 


53 


1999Mar24 


51262.452 


L 


2.0x6.0 


4260-5520 


4.5 


3 


89 


93 


54 


1999Aprll 


51279.533 


LI 


4.2x19.8 


4050-5850 


9 


4 





59 


55 


1999Junl2 


51342.409 


LI 


4.2x19.8 


3999-5796 


9 


3 





29 


56 


1999Junl4 


51344.354 


LI 


4.2x19.8 


4050-5850 


9 


3 


90 


63 


57 


1999Junl6 


51346.335 


LI 


4.2x19.8 


4050-5850 


9 


3 


90 


21 


58 


1999Jull6 


51376.277 


L 


2.0x6.0 


4000-5840 


8 


2 


138 


72 


59 


1999Jul30 


51390.253 


LI 


4.2x19.8 


4000-5800 


9 


2 





34 


60 


1999Aug08 


51399.254 


LI 


4.2x19.8 


4000-5800 


9 


4 





41 


61 


1999Aug09 


51400.245 


LI 


4.2x19.8 


5550-7350 


9 


2 





35 


62 


2000JanlO 


51553.567 


LI 


4.2x19.8 


4000-5800 


9 


3 


90 


78 


63 


2000Jan27 


51570.959 


GH 


2.5x6.0 


4070-7350 


15 


2.5 


90 


83 


64 


2000Jan28 


51571.914 


GH 


2.5x6.0 


4070-7350 


15 


1.5 


90 


72 


65 


2000Febl4 


51588.510 


LI 


4.2x19.8 


4020-5820 


9 


4 


90 


53 


66 


2000Febl5 


51589.506 


LI 


4.2x19.8 


4030-5830 


9 


4 


90 


73 


67 


2000Feb26 


51600.898 


GH 


2.5x6.0 


4560-7850 


15 


2.2 


90 


85 


68 


2000Feb27 


51601.847 


GH 


2.5x6.0 


4310-7600 


15 


2.5 


90 


97 


69 


2000Marl4 


51618.483 


LI 


4.2x19.8 


4050-5850 


9 


4 


90 


34 


70 


2000Marl5 


51619.478 


LI 


4.2x19.8 


4050-5850 


9 


4 


90 


50 


71 


2000Apr01 


51636.414 


LI 


4.2x19.8 


4100-5900 


9 


4 


90 


32 


72 


2000Apr04 


51639.397 


LI 


4.2x19.8 


4050-5850 


9 


3 


90 


41 


73 


2000Apr25 


51659.839 


GH 


2.5x6.0 


4210-7460 


15 


2.2 


90 


61 


74 


2000Apr26 


51660.815 


GH 


2.5x6.0 


4210-7460 


15 


2.2 


90 


48 


75 


2000Mayll 


51676.365 


LI 


4.2x19.8 


4050-5850 


9 


3 


90 


48 


76 


2000Mayl6 


51681.283 


LI 


4.2x19.8 


5550-7350 


9 


3 


90 


52 


77 


2000May20 


51685.321 


LI 


4.2x19.8 


5600-7400 


9 


2 


90 


61 


78 


2000May21 


51686.258 


LI 


4.2x19.8 


4050-5850 


9 


3 


90 


53 


79 


2000May21 


51686.350 


LI 


4.2x19.8 


5550-7350 


9 


2 


90 


41 


80 


2000May25 


51689.786 


GH 


2.5x6.0 


4130-7400 


15 


2.5 


90 


96 


81 


2000May26 


51690.835 


GH 


2.5x6.0 


4134-7380 


15 


2.5 


90 


71 


82 


2000Jun06 


51702.378 


LI 


4.2x19.8 


4020-5820 


9 


3 


90 


56 


83 


2000Jun08 


51704.299 


LI 


4.2x19.8 


4020-5820 


9 


3 


90 


47 


84 


2000Junl5 


51711.308 


LI 


4.2x19.8 


5600-7400 


9 


3 


90 


41 


85 


2000Jun25 


51721.712 


GH 


2.5x6.0 


4690-7980 


15 


2.8 


90 


49 


86 


2000Jull2 


51738.335 


LI 


4.2x19.8 


5560-7360 


9 


3 


90 


33 


87 


2000Jull3 


51739.275 


LI 


4.2x19.8 


5560-7360 


9 


3 


90 


44 


88 


2000Jul28 


51754.269 


LI 


4.2x19.8 


4060-5860 


9 


3 


90 


38 


89 


2000Jul29 


51755.301 


LI 


4.2x19.8 


4040-5840 


9 


2 


90 


46 


90 


2000Jul3031 a 


51756.786 


LI 


4.2x19.8 


4040-5840 


9 


2 


90 


90 


91 


2001Jan25 


51935.482 


LI 


4.2x19.8 


4050-5850 


9 


2 


90 


77 


92 


2001Jan26 


51936.517 


LI 


4.2x19.8 


5550-7350 


9 


2 


90 


58 


93 


2001Jan29 


51938.560 


LI 


4.2x19.8 


4094-5740 


9 


2 


90 


58 
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No 


UT-Date 


JD 


Code 


Aperture 


Sp. range 


Res. 


Seeing 


PA 


S/N 






(2400000+) 




(arcsec) 


(AA) 


(A) 


arcsec 


(deg) 


5160-5220 AA 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


94 


2001Fcb03 


51944.498 


LI 


4.2x19.8 


4094-5740 


9 


2 


90 


54 


95 


2001Febl6 


51957.492 


LI 


4.2x19.8 


4094-5740 


9 


2 


90 


52 


96 


2001Marl3 


51981.867 


GH 


2.5x6.0 


4120-7390 


15 


2.5 


90 


78 


97 


2001Aprl3 


52013.360 


LI 


4.2x19.8 


4094-5740 


9 


3 


90 


69 


98 


2001Aprl3 


52013.451 


LI 


4.2x19.8 


5550-7350 


9 


3 


90 


50 


99 


2001May05 


52034.784 


GH 


2.5x6.0 


4156-7460 


15 


2 


90 


113 


100 


2001Mayl2 


52041.829 


GH 


2.5x6.0 


3600-6900 


15 


1.8 


90 


96 


101 


2001Mayl4 


52043.859 


GH 


2.5x6.0 


3980-7300 


15 


1.8 


90 


125 


102 


2001Junl4 


52074.772 


GH 


2.5x6.0 


4022-7330 


15 


1.8 


90 


93 


103 


2001Junl5 


52075.738 


GH 


2.5x6.0 


4010-7330 


15 


1.8 


90 


100 


104 


2001JullO 


52101.367 


L 


2.0x6.0 


3630-6050 


8 


2 


90 


56 


105 


2001JullO 


52101.367 


L 


2.0x6.0 


5740-8050 


8 


2 


90 


25 


106 


2001Jull2 


52103.283 


L 


2.0x6.0 


3630-6050 


8 


2 





107 


107 


2001Jull2 


52103.283 


L 


2.0x6.0 


5740-8050 


8 


2 





63 


108 


2001Jul21 


52112.340 


LI 


4.2x19.8 


4042-5800 


9 


2 


90 


34 


109 


2001Jul22 


52113.261 


L 


2.0x6.0 


3600-6020 


8 


2 





92 


110 


2001Jul22 


52113.261 


L 


2.0x6.0 


5740-8050 


8 


2 





46 


111 


2001Jul23 


52114.335 


LI 


8.0x19.8 


4094-5746 


9 


2 


90 


47 


112 


2001Jul26 


52117.333 


LI 


8.0x19.8 


4094-5746 


9 


2 


90 


62 


113 


2001Aug09 


52131.300 


LI 


4.2x19.8 


4400-6200 


9 


2 


90 


40 


114 


2002Mayl5 


52410.456 


LI 


4.2x19.8 


4400-6200 


9 


2 


90 


52 


115 


2002Mayl7 


52412.435 


LI 


4.2x19.8 


4400-6200 


9 


2 


90 


49 


116 


2002Jun04 


52430.735 


GH 


2.5x6.0 


3976-7305 


15 


2 


90 


114 



Date 2000 Jul3031 means average of data during two nights; for this case average JD is given. 
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Table B.5. Observed continuum, H^ and H a fluxes. Columns: 1 - UT-Date; 2 - Julian date; 3- a telescope code, 
according to Tabled 4 - F(cont), the continuum flux at 5190 A (in units of 10 -15 ergs -1 cm -2 A -1 ), reduced to the 
1 m telescope aperture 4.2" x 19.8"; 5 - e c , the estimated continuum flux error; 6 - F(Hf3), the H/3 total flux (in units 
of 10~ 13 ergs -1 cm -2 ); 7 - £h/3, the H/3 flux error; 8- F(Ha), the Ha total flux (in units of 10 -13 ergs" 1 cm -2 ); 9 - 
£hq, the Ha flux error. 



UT-Date JD Code F(cont) e F(H /3 ) e F(H„) e 





(2400000+) 




(5190)A 




(4795-5018)AA 




(6500-6800)AA 




1 


2 


3 


4 


5 


6 


7 


8 


9 


1996Janl4 


50097.570 


LI 






7.32 


0.22 


29.90 


1.92 


1996Janl5 


50098.581 


LI 










32.74 


2.10 


1996Febl4 


50127.579 


L 






8.26 


0.24 






1996Febl4 


50127.569 


L 










29.76 


1.98 


1996Febl4 


50128.432 


L 






8.79 


0.25 






1996Febl4 


50128.472 


L 










32.81 


2.18 


1996Marl9 


50162.382 


L 


14.13 


1.03 


7.85 


0.08 






1996Marl9 


50162.399 


L 










28.47 


2.59 


1996Mar21 


50164.395 


L 


12.74 


0.93 


7.89 


0.08 






1996Mar21 


50164.412 


L 










25.04 


2.28 


1996JullO 


50275.283 


L 


16.98 


0.51 


10.04 


0.30 






1996JullO 


50275.353 


L 










41.38 


2.07 


1997Apr05 


50544.414 


LI 


10.83 


0.63 


7.67 


0.38 






1997Apr05 


50544.455 


LI 










36.69 


1.83 


1997Apr08 


50547.331 


L 


11.75 


0.68 


7.14 


0.36 






1997Aprl4 


50553.322 


L 


10.35 


0.31 


7.53 


0.23 






1997Aprl4 


50553.330 


L 










37.53 


.88 


1998Jan21 


50834.632 


L 


14.96 


0.45 


8.96 


0.27 






1998Jan21 


50834.639 


L 










34.65 


1.73 


1998Feb23 


50867.535 


L 


14.36 


0.43 


10.16 


0.30 






1998Feb23 


50867.539 


L 










37.81 


1.89 


1998Apr27 


50931.474 


LI 










42.10 


1.68 


1998Apr30 


50934.438 


LI 










39.81: b 


2.00 


1998Apr30 


50934.467 


LI 


18.49 


0.55 


9.98 


0.30 






1998May04 


50938.289 


L 










36.98 


1.85 


1998May04 


50938.293 


L 


18.85 


0.83 


10.57 


0.12 






1998May06 


50940.433 


L 


20.05 


0.88 


10.4 


0.52 






1998May06 


50940.436 


L 










46.64: 


5.6 


1998May07 


50941.551 


L 


18.92 


0.78 


11.16 


0.28 






1998May08 


50942.452 


L 


18.06 


0.74 


11.72 


0.29 






1998May08 


50942.461 


L 


17.42 


0.71 


11.34 


0.28 






1998May08 


50942.472 


L 










39.60 


1.98 


1998Junl7 


50982.407 


LI 


17.52 


0.18 


11.42 


0.23 






1998Junl9 


50984.287 


LI 


17.26 


0.17 


11.74 


0.23 






1998Junl9 


50984.328 


LI 










44.45 


1.68 


1998Jun20 


50985.281 


LI 










47.31 


1.79 


1998Jun20 


50985.350 


LI 


18.44 


:0.85 


11.55 


0.14 






1998Jun25 


50990.409 


L 


16.41 


0.10 


11.90 


0.17 






1998Jun26 


50991.331 


L 


16.28 


0.10 


11.67 


0.16 






1998Jull4 


51008.778 


GH 


15.25 


0.14 


10.24 


0.19 


42.18 


0.81 


1998Jull7 


51011.683 


GH 


15.05 


0.14 


10.52 


0.20 


43.35 


0.84 


1998Jul26 


51020.703 


GH 


15.25 


0.34 


10.68 


0.31 


41.69 


3.62 
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UT-Datc 


JD 

(2400000+) 


Code 


F(cont) 
(5190)A 


e 


F(Hg) 

(4795-5018)1. 


e 

i 


F(H a ) 
(6500-6800)AA 


e 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1998Jul27 


51021.712 


GH 


15.73 


0.35 


10.26 


0.30 


36.86 


3.20 


1999Janl3 


51191.952 


GH 


13.78 


0.21 


9.41 


0.10 


32.40 


1.04 


1999Janl4 


51193.011 


GH 


14.08 


0.21 


9.34 


0.10 


33.90 


1.08 


1999Jan22 


51200.542 


LI 


13.46 


0.40 


9.44 


0.28 






1999Jan23 


51201.579 


LI 










34.24 


1.10 


1999Jan25 


51203.565 


LI 


14.26 


0.43 


9.62 


0.29 






1999Jan26 


51204.576 


LI 










34.45 


1.10 


1999Marl5 


51252.970 


GH 


15.63 


0.47 


8.57 


0.26 


38.69: 


1.93 


1999Mar24 


51262.444 


LI 


16.59 


0.16 


10.03 


0.16 






1999Mar24 


51262.452 


LI 










36.17: 


1.81 


1999Mar24 


51262.452 


L 


16.53 


0.16 


9.80 


0.16 






1999Aprll 


51279.533 


LI 


15.10 


0.45 


9.66 


0.29 






1999Junl2 


51342.409 


LI 


12.55 


0.25 


9.76 


0.10 






1999Junl4 


51344.354 


LI 


12.66 


0.25 


9.89 


0.10 






1999Junl6 


51346.335 


LI 


13.00 


0.26 


9.81 


0.10 






1999Jull6 


51376.277 


L 


14.06 


0.42 


10.12 


0.30 






1999Jul30 


51390.253 


LI 


10.98 


0.33 


8.36 


0.25 






1999Aug08 


51399.254 


LI 


10.43 


0.31 


8.76 


0.26 






1999Aug09 


51400.245 


LI 










35.52 


1.78 


2000JanlO 


51553.567 


LI 


10.03 


0.30 


9.32 


0.28 






2000Jan27 


51570.959 


GH 


10.47 


0.31 


7.95 


0.26 


33.45 


2.16 


2000Jan28 


51571.914 


GH 


12.52 


0.38 


7.59 


0.25 


36.64 


2.36 


2000Febl4 


51588.510 


LI 


10.55 


0.50 


7.21 


0.17 






2000Febl5 


51589.506 


LI 


11.27 


0.53 


6.97 


0.17 






2000Feb26 


51600.898 


GH 


9.22 


0.13 


6.34 


0.20 


29.62 


2.24 


2000Feb27 


51601.847 


GH 


9.40 


0.13 


6.63 


0.21 


26.60 


2.02 


2000Marl4 


51618.483 


LI 


11.10 


0.33 


7.03 


0.56 






2000Marl5 


51619.478 


LI 


9.98 


0.30 


6.26 


0.50 






2000Apr01 


51636.414 


LI 


8.82 


0.12 


6.60 


0.17 






2000Apr04 


51639.397 


LI 


8.65 


0.12 


6.36 


0.17 






2000Apr25 


51659.839 


GH 


7.83 


0.08 


5.68 


0.06 


28.15 


1.64 


2000Apr26 


51660.815 


GH 


7.73 


0.08 


5.63 


0.06 


25.95 


1.52 


2000Mayll 


51676.365 


LI 


8.59 


0.26 


5.22 


0.16 






2000Mayl6 


51681.283 


LI 










22.84 


1.14 


2000May20 


51685.321 


LI 










21.94 


0.21 


2000May21 


51686.258 


LI 


8.21 


0.25 


5.13 


0.15 






2000May21 


51686.350 


LI 










22.24 


0.21 


2000May25 


51689.786 


GH 


8.13 


0.08 


5.35 


0.18 


22.54 


0.72 


2000May26 


51690.835 


GH 


8.12 


0.08 


5.07 


0.15 


23.59 


1.18 


2000Jun06 


51702.378 


LI 


7.47 


0.22 


5.65 


0.18 






2000Jun08 


51704.299 


LI 


9.31 


0.28 


5.40 


0.17 






2000Junl5 


51711.308 


LI 










19.10 


0.96 


2000Jun25 


51721.712 


GH 


8.09 


0.24 


5.07 


0.15 


18.09 


0.90 


2000Jull2 


51738.335 


LI 










20.26 


0.35 


2000Jull3 


51739.275 


LI 










19.76 


0.35 


2000Jul28 


51754.269 


LI 


9.14 


0.27 


6.37 


0.07 






2000Jul29 


51755.301 


LI 


9.54 


0.29 


6.47 


0.07 






2000Jul3031 a 


51756.786 


LI 


9.01 


0.36 


6.23 


0.17 






2001Jan25 


51935.482 


LI 


8.53 


0.26 


5.85 


0.18 






2001Jan26 


51936.517 


LI 










27.50 


1.38 
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UT-Datc JD Code F(cont) e F(H^) e F(H Q ) e 





(2400000+) 




(5190)A 




(4795-5018)AA 




(6500-6800)AA 




1 


2 


3 


4 


5 


6 


7 


8 


9 


2001Jan29 


51938.560 


LI 


9.58 


0.29 


5.86 


0.18 






2001Feb03 


51944.498 


LI 


8.98 


0.27 


5.90 


0.18 






2001Febl6 


51957.492 


LI 


8.68 


0.26 


6.28 


0.19 






2001Marl3 


51981.867 


GH 


9.09 


0.27 


5.52 


0.17 


23.26 


1.16 


2001Aprl3 


52013.360 


LI 


8.85 


0.27 


4.28 


0.13 






2001Aprl3 


52013.451 


LI 










21.00 


1.05 


2001May05 


52034.784 


GH 


7.71 


0.23 


3.52 


0.11 


17.85 


0.89 


2001Mayl2 


52041.829 


GH 


7.98 


0.08 


3.32 


0.04 


14.97 


0.72 


2001Mayl4 


52043.859 


GH 


7.94 


0.08 


3.34 


0.04 


16.01 


0.77 


2001Junl4 


52074.772 


GH 


9.91 


0.20 


4.55 


0.08 


18.57 


0.82 


2001Junl5 


52075.738 


GH 


9.64 


0.19 


4.44 


0.07 


17.45 


0.77 


2001JullO 


52101.367 


L 


11.27 


0.11 


6.36 


0.22 






2001JullO 


52101.383 


L 










26.01 


0.38 


2001Jull2 


52103.273 


L 










26.56 


0.39 


2001Jull2 


52103.283 


L 


11.31 


0.11 


6.68 


0.23 






2001Jul21 


52112.340 


LI 


11.45 


0.41 


6.27 


0.16 






2001Jul22 


52113.261 


L 


10.65 


0.38 


6.56 


0.17 






2001Jul22 


52113.286 


L 










24.13 


1.21 


2001Jul23 


52114.335 


LI 


11.13 


0.40 


6.28 


0.16 






2001Jul26 


52117.333 


LI 


10.72 


0.30 


6.49 


0.15 






2001Aug09 


52131.300 


LI 


10.54 


0.32 


4.95 


0.15 






2002Mayl5 


52410.456 


LI 


7.32 


0.22 


2.64 


0.08 






2002Mayl7 


52412.435 


LI 


8.11 


0.24 


2.64 


0.08 






2002Jun04 


52430.735 


GH 


7.41 


0.22 


2.42 


0.07 


13.70 


0.68 



a Date 2000 Jul3031 means average of data during two nights; for this case average JD is given. 
b Colon marks unsure value. 



